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Apresentação 


Esse encarte traz dois temas voltados para o controle automático atuante sobre 
o sistema regulador de velocidade e tensão de usinas de energia elétrica onde 
são utilizadas as técnicas do controle moderno no sentido de tornar a operação 
em rede interligada mais confiável. 

O tema principal foca a análise da estabilidade dinâmica do sistema frente a 
ocorrências de abertura de linhas de transmissão quando é verificada a 
comutação entre as condições operativas rede interligada e rede isolada. 

Como é bem sabido, o desempenho de reguladores de velocidade nessas 
situações, dependendo do impacto gerado pela saída de grandes blocos de 
geração, é bastante prejudicado pela alteração significativa no valor do torque 
sincronizante da área restante. 

A proposição de um regulador a base da realimentação de estado é mostrada 
em dois capítulos onde a técnica proposta permite que se ataque esse 
problema de maneira objetiva e prática. 
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CAPÍTULO 1 - ÍNDICES DE DESEMPENHO EM SISTEMAS DE CONTROLE 

AUTOMÁTICO 


RESUMO - Este artigo tem por finalidade apresentar uma discussão acerca da 
utilização de determinados índices de performance e desempenho aplicáveis 
na teoria de controle ótimo de sistemas automáticos. A busca pelo controle 
ótimo, idéia aplicada ao estudo da teoria de controle moderno, baseia-se na 
tomada de decisão com base nesses índices formulados a partir de formas 
quadráticas cujo foco principal é a minimização da energia total gasta pela 
dinâmica desses sistemas no sentido de se obter uma condição de estabilidade 
operativa. A base do estudo levará em conta o princípio de Hamilton a partir 
das equações de Lagrange quando serão implementadas técnicas dentro da 
matemática contínua quanto da matemática discreta. 

Palavras-Chave: Teoria de Controle Ótimo, índices de Desempenho de 
Controle, Controle Moderno. 

Lista das Principais Variáveis 

X : vetor das variáveis de estado de ordem nxl 
U : vetor das variáveis de entrada de ordem mxl 
A : matriz das variáveis de estado da planta de ordem nxn 
B : matriz das variáveis de entrada da planta de ordem nxnn 
F : correspondente da matriz A no plano discreto 
G : correspondente da matriz B no plano discreto 
Q : matriz de ponderação das variáveis de estado de ordem nxn 
Z : matriz de ponderação das variáveis de entrada de ordem mxm 
n : ordem das variáveis de estado do sistema 
m : ordem das variáveis de entrada do sistema 
M : massa 
f :força 

x,y,z,q,p : coordenadas ou variáveis de estado 

1.1. Introdução 

A técnica do controle ótimo aplicada a sistemas de controle automáticos é um 
método ligado ao controle moderno cujo fundamento consiste no controle de 
sistemas pela utilização de uma ferramenta baseada na modelagem por 
variáveis de estado. A abordagem por variáveis de estado tem como 
fundamento a utilização de uma modelagem estruturada na forma matricial com 
grande abrangência nos estudos atuais. 

Assim como em qualquer tratamento de análise de sistemas de controle 
automático, também aqui o importante é a operação estável o que coloca o 
estudo da estabilidade como o fulcro dessa abordagem. 
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Um dos recursos mais importantes tanto na análise quanto na solução do 
problema da instabilidade de sistemas de controle é a utilização da 
realimentação de estado orientada pela minimização de índices de 
desempenho. 

O presente trabalho mostra o desenvolvimento, a elaboração e a aplicação de 
dois tipos de índices, sendo um com foco na dinâmica contínua e outro com 
foco na dinâmica digital dos sistemas de controle. Ver referências [1], [3], [6] e 

m. 

1.2. Revisão da Teoria Clássica 

Qualquer estudo que se faça no campo da teoria de controle envolve com 
conceitos da teoria clássica da dinâmica de partículas e corpos rígidos. Para 
investigar a dinâmica de partículas ou mesmo de corpos rígidos é vantajoso 
especificar as posições do movimento dessas partículas no espaço por meio de 
variáveis independentes do que denotar tais posições através de coordenadas 
retangulares porque nem todas as coordenadas retangulares poderão ser 
independentes durante a análise. Quando isso acontece depara-se com um 
complicador. 

Seja considerar um sistema consistindo de r massas Mi, M2, M3,. dispostas 
no espaço tridimensional cujas coordenadas sejam os seguintes agrupamentos 
(. xiX2,X3); (X 4 ,X 5 ,X 6 ); (x7,xs,xg); ...; (x3r-2,X3r-i,X3r). Seja considerar agora as 
respectivas forças atuantes sobre as r massas como sendo 0 conjunto de 
elementos {UfeM; (kfsje); frfsjg);(Í 3 r- 2 ,f 3 r-i,f 3 r). A equação básica de todo 
movimento será dada pela segunda lei de Newton do seguinte modo: 

Mj Xj = fj i = 1,2,3,... ,3 r (1.1) 

Para que haja concordância nessa proposição, as massas deverão seguir a 
condição: ( Mi=M2=M3 ); ( M 4 =Ms=M 6 ); ( M7=Ms=Mg); ... ; (M 3 r- 2 =M 3 r-i=M 3 r). A 
equação 1.1 representa 0 movimento das massas na ausência de qualquer tipo 
de restrição que possa criar algum tipo de dependência entre as variáveis x,. 
Agora, se 0 sistema admitir as seguintes restrições: 

4 j = 4 j(x v x 2 ,x 3 , ...,x 3r )=o j = 1,2,3,...,s 
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Então, é possível dizer que existem s variáveis linearmente dependentes e que 
o sistema ainda apresenta 3r-s variáveis linearmente independentes de modo 
que: 


4j = 4j(x v x 2 ,x 3 , ...,x 3r )*0 j = (s + 1),(s + 2),(s + 3), ... ,3r 

Assim, para as novas variáveis qi, q 2 , q 3 ,—,qK, onde k=3r-s que formam um 
conjunto de variáveis independentes, tem-se que 

Xj = Xj (q v q 2 , q 3 , ..., q k )*■ ° / = 1,2,3, ... ,3r 


As novas coordenadas qv são chamadas de coordenadas gerais ou 
generalizadas. Para informação, o número de coordenadas gerais é igual ao 
número de graus de liberdade do sistema em questão. Ver referência 3. Para a 
obtenção do princípio de Hamilton, serão utilizadas as equações de Lagrange 
para o movimento em função das coordenadas gerais qi, q 2 , q 3 ,—,qK. 

O princípio de Hamilton é abstraído do princípio de DAIambert que declara que 
a força aplicada a um sistema junto com a força de inércia a ele associada 
formam um conjunto único em equilíbrio, desconsiderando aqui qualquer 
restrição. De maneira geral, pode-se escrever a seguinte expressão 
considerando um espaço físico tridimensional em função do desvio diferencial 
especial <7aqui inserido por comodidade: 


d 2 X: 


I [(Mi 

i = 1 dt 


d 2 yj 


d 2 z ; 


f( x )j) 8 Xj + (Mj ^2 \ y)jí s Vi + (Mj ^2 *( z )i z j] -0 


( 1 . 2 ) 


Segundo o desenvolvimento que pode ser visto na referência [3], a equação 
anterior admite ser realinhada, tomando esta a seguinte forma: 



1 n dX: o dy: o dZ: o O 

- l-gj-) 2 !) + XJ M > s *i + '(y)i Sy i * '(^ Sz i> 


(1.3) 


Na equação 1.3, a partir do segundo membro, notam-se duas parcelas, uma 
representando a energia cinética e a outra a energia potencial. A associação 
considera as n partículas do sistema em questão. A partir daqui não se fará 
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referência aos corpos rígidos por questões de simplificação do texto, nada mais 
que isso. 

Agora é possível sintetizar toda a dinâmica numa única equação em função da 
diferença entre a energia cinética e a energia potencial da seguinte forma: 

n dXj dy: dZ: ,o 

l l/ M í(-à Sx i+ir s y''--á Sz ''>= sí ti( T -vw (i.4) 

Se as posições das partículas do sistema são fixadas para cada tempo 
específico ti e fe, as variações diferenciais 5 x h 5 y, e 5 z, serão nulas o que 
permitirá escrever que: 


s\\ 2 i(T -V)dt = o (1.5) 

A equação 1.5 é conhecida como a equação do princípio de Hamilton e 
determina que a diferença entre a energia cinética Te a energia potencial V de 
um sistema dinâmico é zero na condição de equilíbrio. Na equação 1.5, tem-se 
então que: 


ÕT = 



( 1 . 6 ) 


e 


8 V = -. 1 1 ( f (x)i S x i + f (y)i S Vi + f (z)i sz i) (1 - 7) 

Definindo agora a variável L relativa ao lagrangeano em função do princípio de 
Hamilton tem-se que: 

l = t-v ( 1 . 8 ) 

onde L = ldq j ,q j ) 

Deste modo: 


J/f SL(q j ,q j )dt = 0 (1-9) 

Em termos de diferencial total, a partir da equação 1.9, é possível escrever 
que: 
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( 1 . 10 ) 


ô L(q j ,q j ) = — ôq j + — ôq j i = 1,2,3, ... ,n 
õq, õqi 

Considerando, agora, que o sistema em questão tem n variáveis 
independentes, pela equação 1.9 vem que: 



Ôq j ]dt = 0 
õq, 1 


( 1 . 11 ) 


Sabe-se que a seguinte passagem também é possível: 


8 9 / = 8 


dqj 

dt 


d(8qj) 

dt 


Com isso, a equação 1.11 poderá ser reescrita tomando a seguinte forma: 




( 1 . 12 ) 


Sabe-se que para posições fixas das coordenadas gerais qi, q 2 , q 3 ,—,qn , nos 
tempos específicos ti e fc, as variações diferenciais 5 q, são nulas. Porém, para 
um caso arbitrário em que se possa ter 5 q^O então a seguinte relação deverá 
prevalecer: 


d_ 

dt 


õL . 5L 

õqj J ~dq- 


i = 1,2,3, ... ,n 


(1.13) 


A passagem da equação 1.12 para a equação 1.13 só foi possível quando se 
considerou L variando no domínio do tempo e que as variações das posições 
das novas coordenadas eram diferentes de zero. A equação 1.13 é conhecida 
como a equação de Lagrange para o movimento de partículas dentro do campo 
da dinâmica do ponto material. 

1.3. Equações de Hamilton do Movimento 

As funções de Lagrange podem ser reescritas a partir das equações de 
Hamilton que são dadas por: 


H = T + V 


(1.14) 


Assim 
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L = T -V = 2T - H 


(1.15) 


A proposição agora é trabalhar não só com as variáveis q,, mas também com 
as variáveis p,. As variáveis p, são assim definidas em função das variáveis q,-. 


ÔT 

Pi - — - a 11 q 1 + a 1 2 q 2 + a i3q3 + -+ a i n Qn 
ÕT 

P 2 ~ d^2 ~ 921 Ql + 922 Q2 + 323 q3 + + 32,1 qn 


P3 


ÕT 

õq 3 


~ a 31 Ql + a 32 ^2 + a 33 P3 + 


+ a 3n Qn 


(1.16) 


Pn = 


ÕT 

õQn 


a n1 Ql + a n2 Q2 + a n3 Q3 + — + a nn Qn 


Com a mudança proposta, a energia em jogo será a energia cinética de T que 
representa a energia cinética de todo o sistema de n partículas. Lembrando 
que a energia cinética total, em função da posição de cada partícula do 
sistema, é dada pela seguinte relação: 


7- = ~ S Mjfkjf (1.17) 

2 i = i 

É importante salientar que as variáveis em jogo nesse momento são aquelas 
variáveis linearmente independentes que constituem o sistema uma vez que 
não interessa trabalhar com as variáveis que são linearmente dependentes, 
pois elas não resolvem o problema. 

Aqui é possível associar a equação 1.17 com as variáveis independentes numa 
forma quadrática relacionada à energia cinética envolvida. Assim, 


T 


1 n n n n 

7ÍI ÍI /a-QjQ.; 2T = I [ I 

2 i =1 j = 1 IJ J j = 1 j = 1 


[ajjQjlQi 


(1.18) 


Portanto, com base na equação 1.16, obtém-se: 


2T = 


ÔT 


i = 1 dQj Q ' 


(1.19) 


Aplicando as equações 1.5 e 1.15 relativas ao princípio de Hamilton, a equação 
1.19 poderá ser assim utilizada numa nova formulação, senão veja-se: 
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8 \\ 2 (2T - H)dt = 8 \\ 2 ( Z | f qf/ - H) dt 
r l i = 1 °dj 


( 1 . 20 ) 


Substituindo na equação 1.20 as variáveis p, dadas nas equações 1.16 obtém- 
se a seguinte expressão: 


* (2T - H)dt = 8 f ( z Pi Cfj - H)dt = f 
i = 1 i 


£ [8(p j q j )-5H]dt = 0 
= 1 


( 1 . 21 ) 


Analisando, particularmente, o segundo membro da equação 1.21 vem que: 


8(p j q j ) = p j 8q j + q j 8p j 


e 


u ôH ~ õH „ 

SH = ~8q j Sq i + õp j õp i 


( 1 . 22 ) 


(1.23) 


Substituindo na equação 1.21 as expressões das equações 1.22 e 1.23 e 
agrupando o resultado de maneira adequada obtém-se que: 


tf f (27 - H)dt = f Z [(Pi -~)S Pi + (PjS q; 
h hi = i d Pj 


)8 9 /)] dt = 0 
õqi 


(1.24) 


Trabalhando melhor o segundo membro da equação 1.24, mais precisamente o 
termo \ t t 2 (p j Sq j )dt = o, pode-se assim proceder. Seja a função w = Pjôqj e a sua 

derivada no domínio do tempo conforme o seguinte desenvolvimento: 

Í t (Pi ^ Qi) = Pi^ t (8qi) + 8 q-, £ t (Pj) = P; + ~(Pi)8q- t = Pj 8q,- + p-, 8q,- (1 .25) 

Com base na equação 1.25 pode-se promover a seguinte integração: 

\ t2 (Pi 8 c/j + Pi 8 q s )dt = \ t2 d( p,- 8 q s )= Pj 8 q, f 2 ( 1 ■ 26 ) 

Para os pontos fixos do sistema dinâmico, especificamente, em ti e fe,, 
separadamente, sabe-se que sq j (t 1 ) = sq j (t 2 )=o ■ Deste modo, a partir da 

equação 1.26, obtém-se que: 


(1.27) 


\ t2 (Pi 8 q,-) dt = -\ t2 (Pi 8 qj) dt 
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Considerando-se a expressão obtida da equação 1.27, a equação anterior 1.24 
tomará a seguinte forma: 


r t2 


8\ t2 (2T-H)dt = \\ 2 2 [(q j -^-)õp j -(p j+ ^L)õq j )]dt = 0 
h l ii = 1 B Pj oqj 


(1.28) 


Quando são considerados valores arbitrários para as diferenciais s qj e s qj , tal 
que õqj*o e sqj^o, as seguintes relações ficam evidenciadas: 

i-V.3,-,n (1-29) 

' 'P; 


e 


p:=~™ i =1,2,3,-,n (1.30) 

' õq; 

As equações 1.29 e 1.30 são conhecidas como equações de Hamilton do 
movimento dinâmico para um sistema onde H=T+V, conforme equação 1.14. 

Um sistema com um grau de liberdade /ronde é possível estipular um conjunto 
com um espaço dimensional da ordem de 2k, com as coordenadas gerais q, e 
Pi, i=1,2,3, n , é chamado de espaço de fase. As equações de Hamilton 
fornecem as componentes de velocidade de um determinado ponto do espaço 
de estado. Ficam, então, definidas as equações 1.29 e 1.30 como 
componentes de um conjunto de 2k equações diferenciais de primeira ordem 
na seguinte forma genérica: 

qj=gj(qi,q 2 ’d3. - ,%;Pi,P2’P3 > ->Pk) O- 31 ) 

e 

Pj = hi(di ’Q2’ c l3’ ■■■’%'> PvP2’P3’ - ’Pk) O- 32 ) 

i = 1,2,3, -,/c 

Desde que essas 2k variáveis sejam a menor coleção de números que possam 
ser especificadas no sentido de descrever completamente a dinâmica de um 
determinado sistema essas variáveis são as variáveis de estado desse 
sistema. 
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1.4. Aplicação na Teoria de Controle 


A teoria de controle moderno se utiliza dos chamados índices de desempenho 
para a efetivação de sua ação sobre a dinâmica dos sistemas em geral. Esses 
índices nada mais são que formas quadráticas em função da combinação das 
diversas variáveis de estado que compõem o sistema. Elas formam um 
conjunto mínimo de equações que descrevem completamente a dinâmica do 
sistema em questão. 

A aplicação do princípio de Hamilton da teoria clássica da dinâmica, assim 
como foi visto anteriormente, introduz o conceito de espaço de fase estendido 
ao conceito de espaço de estado com um forte apelo na formulação do 
conceito de controle ótimo com base em índices desempenho. 

A diferença é que não só serão utilizadas as variáveis de estado, mas também 
as variáveis de entrada externa ao sistema. 

Existem várias proposições de formulação dos índices de desempenho, 
também conhecidos como índices de performance. 

Um exemplo bem usual é aquele que agrega as variáveis de estado e ao 
mesmo tempo as variáveis de entrada externa na forma mostrada a seguir, 
onde T denota matriz transposta: 

l= 1 -\ t 0 [X T QX + U T ZU]dt (1-33) 

relativamente à equação de uma planta ou sistema de controle automático na 
conhecida forma canônica de variáveis de estado, representada em seguida. 

X = AX + BU (1-34) 

A técnica de realimentação de estado, um procedimento comum no processo 
de estabilização em sistemas de controle automático, pode ser equacionada 
dentro da forma padrão de variáveis de estado a partir da equação matricial da 
planta 1.34. Pala realimentação de estado tem-se que: 


u = -Kr x 


(1.35) 
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onde K r é a matriz de realimentação de ordem qxn. Deste modo, considerando 
a equação 1.34, obtém-se: 


X = [ A-BKr] X (1.36) 

Em termos do índice de desempenho dado na equação 1.33, a equação 
tomará o seguinte aspecto: 

l = ^\ t 0 [X T QX + (-K r X) J Z(-K r X)]dt = ^\ t 0 X T [Q + K r T ZK r ] Xdt C 1 -37) 

Pode-se facilmente concluir que o índice de desempenho, aqui formulado, 
agrupa três importantes características que são básicas na pesquisa para a 
obtenção de uma matriz de realimentação de estado, a lei de transferência de 
um estado inicial para um estado final, as energias externas aplicadas ao 
sistema que garantem essa transferência e o período de tempo em que são 
processados os eventos. O sistema deverá ser estável nesse período. 

Para o presente problema, a função do hamiltoniano conterá as variáveis 
independentes x, que serão as próprias variáveis de estado do sistema além de 
contar com as variáveis p, de ponderação, para i=1,2,3, ... , n; onde n é a 
ordem máxima do sistema. De maneira compacta, o hamiltoniano terá a 
seguinte forma: 


H = H(x j ,u k ,x j ) 


(1.38) 


com ' - 12,3, ■■■ ,n 

k = 1,2,3, -,m 

Os textos clássicos, referências [1], [3] e [4] sugerem que o hamiltoniano, 
considerando o conjunto de equações na forma canônica de variáveis de 
estado, tenha a seguinte formulação: 

H = ^X T QX + ^U T ZU + P T X = ^X T QX + ^U T ZU + P T [AX + BU] 0 -39) 

Considerando o princípio de Hamilton, conforme descrito pelas equações 1.29 
e 1.30, é possível, tendo em vista a equação 1.39, a obtenção das seguintes 
expressões: 
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(1.40) 


— = — P T [AX + BU] = AX + BU=X 
ÕPõP 


e 


-™ = -±d X T QX] _± [p T AX] = p (1.41) 

ÕX ÔX l 2 j ÕX l 1 

Será possível, também, a utilização de uma terceira equação diferencial 
envolvendo as variáveis de entrada. 


ÔH 

ÔU 


— U T ZU] + — [P T BU] = 0 

ÕU 2 J ÕU 


(1.42) 


Pelo estudo em que se procura a definição de um índice de desempenho mais 
realista possível, como aquele fornecido pela equação 1.33, a introdução dos 
critérios de estabilidade de sistemas de controle segundo Lyapunov, se faz 
necessária e sua aplicação pode ser explicada a seguir. Ver referências [1] e 
[ 6 ], 

Derivando a equação 1.33, relativa ao índice de desempenho, no domínio do 
tempo, tem-se que: 


V(X)=^[X T QX + U T ZU] (1-43) 

Com a inclusão dos efeitos das próprias variáveis de estado no cômputo do 
cálculo variacional é possível escrever a função matricial F dando a seguinte 
relação: 

F=V(X)+P T X (1-44) 

Comparando as equações 1.39 e 1.44, observa-se que a função Fnada mais é 
que o hamiltoniano onde F=H. 

Convém lembrar que a função hamiltoniano de um sistema dinâmico está 
associada ao balanço entre as energias cinética e potencial e que dentro do 
escopo da teoria de controle, no plano de estado, representa o balanço geral 
da energia intrínseca desse sistema. 
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A utilização de um índice de desempenho na forma do hamiltoniano visa 
trabalhar com a ideia de controle sob a óptica de níveis mínimos de energia. 
Nesse sentido, essa ideia também evidencia a questão da operação estável. 


A equação 1.42 permite que se promova o seguinte desenvolvimento: 


- = l-U T ZU] + — [P T BU] = u T z+ p t b=o 

ôu eu 2 õu 


(1.45) 


(U T Z) T + (P T B) T = 0 Z T U = -B T P U = -ZT 1 B T P 


onde adota-se Z ? =Z 


Comumente na literatura de controle automático, a matriz P é uma matriz 
simétrica que varia no domínio do tempo contendo a seguinte forma: 

p = kx (1-46) 

Deste modo, o vetor de variáveis de entrada da equação 1.45 tomará o 
seguinte aspecto: 

U = -ZT 1 B r KX :.U = rX (1.47) 


onde r = -zr 1 B r K 

Derivando a equação 1.46 em relação ao tempo, e já considerando a matriz K 
uma matriz de elementos constantes, o que é mais comum, tem-se: 

p = kx (1.48) 

Promovendo o devido desenvolvimento da equação 1.41 anterior, tem-se a 
seguinte expressão: 

-X T Q-P T A-P = 0 (X T Q) T +(P T A) T + (P) T = 0 

Levando-se em conta a condição de que (p) T = p, e que a matriz Q é simétrica, 
para a equação anterior vem que: 

(X T Q) T + (P T A) T + P = 0 Q T X+A T P + P = 0 


logo 
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ou melhor 


q t x+a t p+p = qx+a t p+kX = o 


QX+A T P+KX = QX + A T P+K[AX+BU] = 0 


ou então 


KAX + KBU + QX + A 7 P=0 (1.49) 

Trazendo-se na equação 1.49 o vetor das variáveis U da equação 1.47, 
finalmente obtém-se a conhecida equação de Riccati que permite a 
consecução da matriz de reaiimentação K do sistema sob a perspectiva do 
controle ótimo. Assim sendo, tem-se que: 

kax + a t kx+qx-kbz~ 1 b t kx = o 


finalmente 


KA + A r K+Q-KBZr 1 B J K=0 (1.50) 

Para o caso em que a matriz K seja composta por funções dependentes do 
tempo, genericamente, ter-se-ia que: 

k+KA + A T K+Q-KBZr 1 B T K=0 (1.51) 

Como observação, vale lembrar que a solução das equações 1.50 e 1.51 
deverá ser obtida por intermédio de recursos computacionais, pois não se trata 
de operação analítica. 

1.5. Abordagem Discreta 

Assim como foi visto anteriormente dentro do campo da matemática contínua, o 
processo de discretização ou discreto também admite que se utilize da forma 
canônica de variáveis de estado, apenas que agora a abordagem é adaptada 
aos processos onde o controle é promovido por computador. As variáveis 
agora se tornam números para serem processados através de fórmulas 
puramente matemáticas ao invés de grandezas físicas tangíveis. A forma 
discreta da equação 1.34 pode ser assim descrita então: 
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X(kz + 1) = F(kz)X(kz)+G(kz)U(kz) 


(1.52) 


COm k = 1,2,3, ...n 

onde ré o passo de integração adotado no processo de discretização. Ver 
referência [4], [5] e [8]. 

O índice de desempenho aqui proposto tem semelhança com o índice 
formulado para o caso em que o modelo estava no plano da matemática 
contínua e pode ser assim descrito: 

IP= ÍH(k) (1.53) 

k=0 

Assim, através do processo de minimização do índice de desempenho discreto 
dado pela equação 1.53, é possível implementar a ação de otimização 
operativa para o sistema contínuo que foi discretizado. 

As formulações que reúnem as condições básicas a serem obedecidas para 
que sejam garantidos os objetivos de minimização do índice de desempenho 
devem seguir as equações matriciais seguintes, bem semelhantes àquelas que 
foram adotadas para a análise no plano contínuo. 

Seguindo o princípio de Hamilton, tem-se que: 

= o=FX(k)+GU(k)-X(k+1) (1-54) 

õ P 


e 


^!íhl = p(k) = QX(k)+F T P(k + 1) (1.55) 

õX 

Será possível, também, a utilização de uma terceira equação diferencial 
envolvendo as variáveis de entrada. 

^t!íhl = o = ZU(k)+ G 7 " P(k+i) (1-56) 


com k = 0 , 1 , 2 , ...(n- 1 ) 
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A ideia de utilização dos índices ditados pelas equações 1.54, 1.55 e 1.56, 
ligeiramente modificados em relação àqueles propostos nas equações 1.40, 
1.41 e 1.42 na abordagem contínua, não foge muito da ideia de minimização. 
Ver referências [1] e [4], 

As seguintes relações podem ser explicitadas a partir daqui: 

X(k + 1) = FX(k) + GU(k) (1-57) 


e 


P(k) = Q X(k)+F T P(k + 1) 


(1.58) 


e também 


U(k) = -Z~ 1 G T P(k + 1) 


(1.59) 


com k = 0,1,2, ...(n- 1 ) 

Substituindo a expressão de U(k) da equação 1.59 na equação 1.57 vem que: 

X(k + 1)=FX(k)-GZ~ 1 G T P(k + 1) (1.60) 

A matriz P, também conhecida como multiplicador de Lagrange, tem uma 
relação com o vetor das variáveis de estado do sistema numa relação 
semelhante à relação da equação 1.46 quando da abordagem contínua do 
problema. Assim, pode-se escrever numa abordagem discreta que: 

P(k) = K(k)X(k) (1.61) 

Igualando as relações das equações 1.58 e 1.61 obtém-se que: 

K(k)X(k) = QX(k) +F T P(k + 1 ) = QX(k) +F T K(k + 1)X(k + 1 ) (1-62) 

Da equação 1.60 obtém-se que: 

X(k + 1)=FX(k)-GZ~ 1 G T K(k + 1)X(k + 1) (1.63) 
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Da equação 1.63, anterior, isola-se a expressão para a matriz X(k+1). Assim, 
tem-se que: 

X(k + 1) = [l + GZ~ 1 G T Kfk + I)]- 1 FX(k) (1-64) 

onde lé a matriz identidade 

Levando-se a expressão de X(k+1) da equação 1.64 na equação 1.62 vem 
que: 

K(k)X(k ) = QX(k) + F T K(k + 1 )[l + GZ~ 1 G T K(k + 1)]~ 1 FX(k) 

Eliminando-se o vetor variáveis de estado X entre os membros da equação 
precedente vem: 


K(k) = Q+ F T K(k + 1)[l + GZ~ 1 G T K(k + 1)]~ 1 F (1.65) 

A equação 1.65 é a equação de Riccati em sua forma discreta que realmente 
contém a relação iterativa possibilitando, através de uma implementação 
computacional, que se obtenha a matriz K. 

De posse da matriz K, através de um simples procedimento matemático, 
obtém-se a matriz P(k+1) que, levada na equação 1.59, fornecerá a matriz das 
variáveis de realimentação da equação U(k). 

Da equação 1.58 extrai-se a expressão da matriz P(k+1) da seguinte maneira: 

P(k + 1) = [F T ]~ 1 [P(k)-QX(k)] (1.66) 

Levando-se a expressão de P(k+1) da equação 1.66 na equação 1.59 vem 
que: 


U(k ) = -Z~ 1 G T [F T r 1 [P(k)-Q X(k)] 

Levando, agora, na equação anterior, a expressão de P(k) da equação 1.61, 
obtém se: 

U(k) = -Z~ 1 G T [F T r UK(k)X(k)-QX(k)] = ~Z~ 1 G J [F T ]~ 1 [K(k)-Q]X(k) = í X(k) (1-67) 
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A equação 1.67 reúne, a partir da determinação da matriz de realimentação F, 
os elementos que garantem o controle ótimo via realimentação de estado. 

Para voltar ao plano contínuo bastará utilizar a fórmula seguinte: 

r (t)=r(k) ( 1 . 68 ) 


1.6. Considerações Finais 

O procedimento iterativo utilizado para a computação digital da matriz K da 
equação de fíiccati terá seu sucesso garantido a partir do momento em que 
houver convergência dos valores sempre a partir de condições iniciais 
favoráveis. 

Algumas observações se fazem necessárias: 

a) Durante o processo de cálculo da matriz da planta no plano discreto, 
matriz F, deve-se assegurar que o passo de integração adotado seja no 
máximo 10% inferior a menor constante de tempo envolvida no processo. 

b) Cuidados especiais deverão ser tomados com subsistemas com laços de 
realimentação. 

c) A condição de estabilidade aqui adotada consiste em considerar, a partir 
do prévio conhecimento das equações fásicas na forma canônica de 
variáveis de estado, que um determinado sistema somente será declarado 
estável se, e somente se, todos os seus autovalores associados à matriz da 
planta do modelo contínuo não apresentarem partes reais positivas. 

d) A matriz de ponderação relativa às variáveis de entrada, matriz Z no plano 
discreto, se tomada como diagonal unitária será suficiente para o bom 
desenvolvimento do processo computacional. 

e) A matriz Q de ponderação relativa às variáveis de estado tem importância 
capital no desempenho do processo. Devem-se testar várias matrizes a 
começar por uma matriz diagonal, porém escolher e testar os valores antes. 

f) A matriz K obtida a partir da resolução da equação de Riccati será sempre 
simétrica desde que a matriz de ponderação Q o seja. 
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g) Partir sempre da condição em que a matriz Q de ponderação seja uma 
matriz definida positiva, no mínimo semi-definida positiva. 
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CAPÍTULO 2 - REGULADOR DE VELOCIDADE POR CONTROLE 
MODERNO VIA EQUAÇÃO DE RICCATI 


RESUMO - Este artigo tem por finalidade apresentar o desenvolvimento e a 
aplicação de um controle suplementar ao regulador de velocidade de uma 
turbina hidráulica por ação de um controlador do tipo discreto baseado em 
minimização de índices de desempenho com base no controle ótimo, uma 
técnica do controle moderno com a abordagem do problema por variáveis de 
estado. Será utilizado o método de realimentação de estado aplicado a uma 
situação de instabilidade sobre o controlador PID tradicional. O sistema de 
controle e o processo a ser utilizado pertencem a uma instalação real cujo 
modelo foi validado previamente através de ensaios de campo. 

Palavras-Chave: Teoria de Controle Ótimo, Reguladores de Velocidade, 
Controle Moderno. 

Lista das Principais Variáveis 

X : vetor das variáveis de estado [nxl] 

U : vetor das variáveis de entrada [mxl] 

A : matriz das variáveis de estado da planta [nxn] 

B : matriz das variáveis de entrada da planta [nxm] 

F : correspondente da matriz A no plano discreto[nxn] 

G : correspondente da matriz B no plano discreto[nxm] 

Q : matriz de ponderação das variáveis [nxn] 

Z: matriz de ponderação das variáveis de entrada [nxn] 

F : matriz de realimentação de estado [nxl] 
n : ordem das variáveis de estado do sistema 
m : ordem das variáveis de entrada do sistema 
M : massa 

A Pm(S) : variação da potência mecânica em pu 

APe(S) : variação da potência mecânica em pu 

A co(S) : variação da velocidade em pu 

A ô(S) : variação do ângulo elétrico de potência em radianos 

A Y(S) : variação do curso do servomotor principal em pu 

A z(S) : variação do spool da válvula piloto em pu 

U(1) : entrada externa de sinais de realimentação 

X(1) : variável de controle integrador em pu 

X(2) : curso da válvula piloto em pu 

X(3) : curso do servomotor principal em pu 

X(4) : variável associada à tubulação forçada em pu 

X(5) : velocidade ou freqüência em pu 

X(6) : ângulo elétrico de potência em radianos 

bp : estatismo permanente em pu/pu 

Td : constante de tempo do estatismo transitório em s 

Td : constante de tempo de washout do regulador em s 

KLP : constante de proporcionalidade sem dimensão 

Ti : constante de tempo associada à válvula proporcional em s 

Ty : constante de tempo do servomotor principal em s 

Ymáx : curso máximo do servomotor em pu 

Ymín : curso mínimo do servomotor em pu 
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Tw : constante de inércia da água em s 
2H : constante de inércia do conjunto girante em s 
D : coeficiente de amortecimento da carga em pu/pu 
T : coeficiente de torque sincronizante em pu/radianos 
Cl : comutador de realimentação de estado 
C2 : comutador de rede on-line para rede off-line 

2.1. Introdução 

A técnica de controle ótimo aplicada a sistemas de controle automático é um 
método ligado ao controle moderno cuja base consiste no controle de sistemas 
pelo uso de uma ferramenta baseada na modelagem por variáveis de estado. A 
abordagem por variáveis de estado é focada no uso de um modelo estruturado 
na forma matricial. 

Como em qualquer tratamento de análise de sistemas de controle automático, 
também aqui o importante é a operação estável fato que coloca o estudo da 
estabilidade como o ponto fundamental dessa abordagem. 

Uma das características mais importantes na análise e solução do problema de 
instabilidade do sistema de controle é o uso da realimentaçaão de estado 
guiado pela minimização dos índices de desempenho. O presente trabalho 
mostra a aplicação desse método ao modelo, já devidamente validado, de um 
sistema regulador de velocidade de uma turbina hidráulica da UHE 
Emborcação (4 x 298 MW) localizada no Triângulo Mineiro. Ver referência [1]. 

2.2. Revisão da Teoria Clássica 

Qualquer estudo que se faça no campo da teoria de controle envolve com 
conceitos da teoria clássica da dinâmica de partículas e corpos rígidos. Para 
investigar a dinâmica de partículas ou mesmo de corpos rígidos é vantajoso 
especificar as posições do movimento dessas partículas no espaço por meio de 
variáveis independentes do que denotar tais posições através de coordenadas 
retangulares porque nem todas as coordenadas retangulares poderão ser 
independentes durante a análise. Quando isso acontece depara-se com um 
complicador. 

O princípio de Hamilton é abstraído do princípio de D'Alambert que declara que 
a força aplicada a um sistema junto com a força de inércia a ele associada 


24 



formam um conjunto único em equilíbrio, desconsiderando aqui qualquer 
restrição. 


sf 1 (T-V)dt=0 ( 2 . 1 ) 

A equação 2.1 é conhecida como a equação do princípio de Hamilton e 
determina que a diferença entre a energia cinética Te a energia potencial V de 
um sistema dinâmico é zero na condição de equilíbrio. 

Definindo agora a variável L relativa ao lagrangeano em função do princípio de 
Hamilton tem-se que: 


L = T -V 


( 2 . 2 ) 


onde L = L(q j ,q j ) 


Deste modo: 


\ t t 2 1 SL(q j ,q j )dt = 0 (2.3) 

Em termos de diferencial total, a partir da equação 2.3, é possível escrever 
que: 


SL(q„q,) = ^Sq, + ^Sq, i = 1,2,3 ,... ,n (2.4) 

õQj dqj 

Considerando, agora, que o sistema em questão tem n variáveis 
independentes, pela equação 2.3 vem que: 


Jf? I [—ôq j + — ôq j ]dt = o 

!t1 i = 1 õ qj cq j 


(2.5) 


Com isso, a equação 2.5 poderá ser reescrita tomando a seguinte forma: 


r t2, " r d ôL dL 

y,l/di < ãí; >+ Wi ISq i )d,=0 


( 2 . 6 ) 


Sabe-se que para posições fixas das coordenadas gerais qi, q 2 , q 3 ,—,qn, nos 
tempos específicos ti e t 2 , as variações diferenciais ôqi são nulas. Porém, para 
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um caso arbitrário em que se possa ter 5 q^O então a seguinte relação deverá 
prevalecer: 


d_ ÕL_ _ õL 

dt õc/j ' Õqi 


i = 1,2,3, ...,n 


(2.7) 


A passagem da equação 2.6 para a equação 2.7 só foi possível quando se 
considerou L variando no domínio do tempo e que as variações das posições 
das novas coordenadas eram diferentes de zero. A equação 2.7 é conhecida 
como a equação de Lagrange para o movimento de partículas dentro do campo 
da dinâmica do ponto material. 

2.3. Equações de Hamilton do Movimento 

As funções de Lagrange podem ser reescritas a partir das equações de 
Hamilton que são dadas por: 


H=T + V 


( 2 . 8 ) 


Assim 


L = T -V = 2T — H 

A proposição agora é trabalhar não só com as variáveis q h mas também com 
as variáveis p,-. As variáveis p, são assim definidas em função das variáveis c/j . 


õT 


Pi - — - a 11 q 1 + a 12 q 2 + a 13 q 3 + ... 

+ a 1n dn 

ÕT 


P2 ~ ~ a 2ldl + a 22 d2 + a 23 % + 

... +a 2n dn 

ÕT 


P3 - ~ ~ a 3ldl + a 32 d2 + a 33 % + 

... +a 3n q n 

ÕT 


Pn ~ ~ a n1 dl + a n2 d2 + a n3 d3 + 

... +a nn q n 


Com a mudança proposta, a energia em jogo será a energia cinética de T que 
representa a energia cinética de todo o sistema de n partículas. Lembrando 
que a energia cinética total, em função da posição de cada partícula de massa 
M do sistema, é dada pela seguinte relação: 
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( 2 . 10 ) 


1 n 2 
T = ~ X Mj (Xjj 

2 i = 1 

É importante salientar que as variáveis em jogo nesse momento são aquelas 
variáveis linearmente independentes que constituem o sistema uma vez que 
não interessa trabalhar com as variáveis que são linearmente dependentes, 
pois elas não resolvem o problema. Aqui é possível associar a equação 2.10 
com as variáveis independentes numa forma quadrática relacionada à energia 
cinética envolvida. Assim: 

T = X [ X [ajjàijjcij} 2T = X [ X [aj:q:]qj (2.11) 
* i = 1 j = 1 11 i — 1 j = 1 11 

Portanto, com base na equação 2.9, obtém-se: 

n õT 

2T = ,l^ (2 ' 12) 

Aplicando as equações 2.1 e 2.8 relativas ao princípio de Hamilton, a equação 
012 poderá ser assim utilizada numa nova formulação, senão veja-se: 

ô\ t2 (2T-H)dt = ô\ t2 ( X qj-H)dt (2.13) 

7 1 i = 1 tíC li 

Substituindo na equação 2.13 as variáveis p,dadas nas equações 2.9 obtém-se 
a seguinte expressão: 

8íj*(2T - H)dt = ,) \ t2 ( i Pi q. -H)dt = \\ 2 X ^ [S( Pj q^-õH]dt = 0 (2.14) 

Analisando, particularmente, o segundo membro da equação 2.14, uma análise 
mais depurada que pode ser vista na referência [3] conduz a seguinte relação: 

8 ff 2 (2T - H)dt = J.‘ 2 X [(q j - — )8p j -(p i + - jSq,)] dt = 0 (2.15) 

7 7 i = 1 °Pj õqj 

Quando são considerados valores arbitrários para as diferenciais 8 q, e 8 q ,, 
tal que 8q, * 0 e 8 q, * 0 , as seguintes relações ficam evidenciadas: 
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ÕH 

dpi 


i = 1,2,3, —,n 


(2.16) 


Qi = 


e 


Pi=~— i = 1,2,3, ••• ,n (2.17) 

õq-, 

As equações 2.16 e 2.17 são conhecidas como equações de Hamilton do 
movimento dinâmico para um sistema onde H=T+V, conforme equação 2.7. 

Um sistema com um grau de liberdade /conde é possível estipular um conjunto 
com um espaço dimensional da ordem de 2k, com as coordenadas gerais g, e 
Pi, i=1,2,3,n , é chamado de espaço de fase. 

As equações de Hamilton fornecem as componentes de velocidade de um 
determinado ponto do espaço de estado. Ficam, então, definidas as equações 
2.16 e 2.17 como componentes de um conjunto de 2k equações diferenciais de 
primeira ordem na seguinte forma genérica: 

Qi = g i (q 1 ,q 2 ,q 3 , - ,q k ! P 1 ,P 2 ,P 3 ,->Pk) ( 2 - 18 ) 

e 

Pi = h i (q 1 ,q 2 ,q 3i ...,q k ;p 1 ,p 2 ,p 3 ,...,p k ) (2-19) 

/ = 1,2,3, -,/c 

Desde que essas 2k variáveis sejam a menor coleção de números que possam 
ser especificadas no sentido de descrever completamente a dinâmica de um 
determinado sistema essas variáveis são as variáveis de estado desse 
sistema. 

2.4. Aplicação na Teoria de Controle 

A teoria de controle moderno se utiliza dos chamados índices de desempenho 
para a efetivação de sua ação sobre a dinâmica dos sistemas em geral. Esses 
índices nada mais são que formas quadráticas em função da combinação das 
diversas variáveis de estado que compõem o sistema. Elas formam um 
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conjunto mínimo de equações que descrevem compietamente a dinâmica do 
sistema em questão. 

A aplicação do princípio de Hamilton da teoria clássica da dinâmica, assim 
como foi visto anteriormente, introduz o conceito de espaço de fase estendido 
ao conceito de espaço de estado com um forte apelo na formulação do 
conceito de controle ótimo com base em índices desempenho. 

A diferença é que não só serão utilizadas as variáveis de estado, mas também 
as variáveis de entrada externa ao sistema. 

Existem várias proposições de formulação dos índices de desempenho, 
também conhecidos como índices de performance. 

Um exemplo bem usual é aquele que agrega as variáveis de estado e ao 
mesmo tempo as variáveis de entrada externa na forma mostrada a seguir, 
onde T denota matriz transposta: 

I = ^\ t 0 [X T QX + U T ZU]dt (2.20) 

relativamente à equação de uma planta ou sistema de controle automático na 
conhecida forma canônica de variáveis de estado, representada em seguida. 

X = AX + BU (2.21) 

A técnica de realimentação de estado, um procedimento comum no processo 
de estabilização em sistemas de controle automático, pode ser equacionada 
dentro da forma padrão de variáveis de estado a partir da equação matricial da 
planta 2.21. 

Pala realimentação de estado tem-se que: 


u = -Kr x 


( 2 . 22 ) 


onde K r é a matriz de realimentação de ordem qxn. 

Deste modo, considerando a equação 2.21, obtém-se: 

X = [A-BKr] X (2.23) 
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Em termos do índice de desempenho dado na equação 2.20, a equação 
tomará o seguinte aspecto: 

l = ^\o [xTQX+ (- K r X ) Tz (~ K r X)]dt= 1 -\ t 0 X T [Q + K r T ZK r ]Xdt (2.24) 

Pode-se facilmente concluir que o índice de desempenho, aqui formulado, 
agrupa três importantes características que são básicas na pesquisa para a 
obtenção de uma matriz de realimentação de estado, a lei de transferência de 
um estado inicial para um estado final, as energias externas aplicadas ao 
sistema que garantem essa transferência e o período de tempo em que são 
processados os eventos. O sistema deverá ser estável nesse período. 

Para o presente problema, a função do hamiltoniano conterá as variáveis 
independentes x, que serão as próprias variáveis de estado do sistema além de 
contar com as variáveis p, de ponderação, para i=1,2,3, ... , n; onde n é a 
ordem máxima do sistema. 

De maneira compacta, o hamiltoniano terá a seguinte forma: 

H = H(x j ,u k ,x j ) (2.25) 


com 


/ = 1,2,3, ■ ■■,n 
k = 1,2,3, ■■■,m 


Os textos clássicos, referências [1], [3] e [4] sugerem que o hamiltoniano, 
considerando o conjunto de equações na forma canônica de variáveis de 
estado, tenha a seguinte formulação: 


H = - X T QX + -u T ZU + P T X = 

2 2 

= 1 - 2 X T QX+ 1 -U T ZU + P T [AX + BU] 


(2.26) 


Considerando o princípio de Hamilton, conforme descrito pelas equações 2.16 
e 2.17, é possível, tendo em vista a equação 2.26, a obtenção das seguintes 
expressões: 


— = — P T [AX+BU] = AX+BU=X (2.27) 

õ P õ P 
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e 


-Tx-TX ÍL 2 XTQXI -TX ÍPTAX1 - P < 2 ' 28) 

Será possível, também, a utilização de uma terceira equação diferencial 
envolvendo as variáveis de entrada. 

— = —[ 1 - U T ZU) + —[P T BU] = 0 (2.29) 

õU õU 2 qu 

Pelo estudo em que se procura a definição de um índice de desempenho mais 
realista possível, como aquele fornecido pela equação 2.20, a introdução dos 
critérios de estabilidade de sistemas de controle segundo Lyapunov, se faz 
necessária e sua aplicação pode ser explicada a seguir. Ver referências [1] e 
[ 6 ], 

Derivando a equação 2.20, relativa ao índice de desempenho, no domínio do 
tempo, tem-se que: 


V(X)=^[X t QX + U t ZU] (2.30) 

Com a inclusão dos efeitos das próprias variáveis de estado no cômputo do 
cálculo variacional é possível escrever a função matricial F dando a seguinte 
relação: 


F =V(X)+ P T X (2.31) 

Comparando as equações 2.26 e 2.31, observa-se que a função F nada mais 
é que o hamiltoniano onde F=H. 

Convém lembrar que a função hamiltoniano de um sistema dinâmico está 
associada ao balanço entre as energias cinética e potencial e que dentro do 
escopo da teoria de controle, no plano de estado, representa o balanço geral 
da energia intrínseca desse sistema. 

A utilização de um índice de desempenho na forma do hamiltoniano visa 
trabalhar com a ideia de controle sob a óptica de níveis mínimos de energia. 
Nesse sentido, essa ideia também evidencia a questão da operação estável. 
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2.5. Abordagem Discreta 


Assim como foi visto anteriormente dentro do campo da matemática contínua, o 
processo de discretização ou discreto também admite que se utilize da forma 
canônica de variáveis de estado, apenas que agora a abordagem é adaptada 
aos processos onde o controle é promovido por computador. As variáveis 
agora se tornam números para serem processados através de fórmulas 
puramente matemáticas ao invés de grandezas físicas tangíveis. 

A forma discreta da equação 2.32 pode ser assim descrita então: 

X(kT + 1) = F(kT)X(kT) + G(kT)U(kT) (2.32) 

com k = 1,2,3, ...n 

onde ré o passo de integração adotado no processo de discretização. Ver 
referência [4], [5] e [8]. 

O índice de desempenho aqui proposto tem semelhança com o índice 
formulado para o caso em que o modelo estava no plano da matemática 
contínua e pode ser assim descrito: 

P = (2-33) 

Assim, através do processo de minimização do índice de desempenho discreto 
dado pela equação 2.33, é possível implementar a ação de otimização 
operativa para o sistema contínuo que foi discretizado. 

As formulações que reúnem as condições básicas a serem obedecidas para 
que sejam garantidos os objetivos de minimização do índice de desempenho 
devem seguir as equações matriciais seguintes, bem semelhantes àquelas que 
foram adotadas para a análise no plano contínuo. 

Seguindo o princípio de Hamilton, tem-se que: 

^fíhl = 0 =FX(k)+GU(k)~ X(k + 1) (2.34) 


e 
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(2.35) 


dH(k) t 

—L -t- = P(k) = QX(k) + F' P(k + 1) 
õX 

Será possível, também, a utilização de uma terceira equação diferencial 
envolvendo as variáveis de entrada. 


= o = ZU(k)+ G T P(k + 1 ) (2.36) 

ÕU 


com k = 0,1,2, ...(n-1) 

A ideia de utilização dos índices ditados pelas equações 2.34, 2.35 e 2.36, 
ligeiramente modificados em relação àqueles propostos nas equações 2.27, 
2.28 e 2.29 na abordagem contínua, não foge muito da ideia de minimização. 
Ver referências [1] e [4], 

As seguintes relações podem ser explicitadas a partir daqui: 

X(k+1) = FX(k) + GU(k) (2.37) 


e 


P(k) = Q X(k)+F T P(k + 1) 


(2.38) 


e também 


U(k) = -z 1 G r P(k + i) 


(2.39) 


com k = 0,1,2, ...(n-1) 

Substituindo a expressão de U(k) da equação 2.39 na equação 2.37 vem que: 

X(k + 1)= FX(k)-Gz~ 1 G r P(k + i) (2.40) 

A matriz P, também conhecida como multiplicador de Lagrange, tem uma 
relação com o vetor das variáveis de estado do sistema numa relação 
semelhante à relação da equação quando da abordagem contínua do 
problema. Assim, pode-se escrever numa abordagem discreta que: 

P(k) = K(k)X(k) (2.41) 
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Igualando as relações das equações 2.38 e 2.41 obtém-se que: 

K(k)X(k) = OX(k)+ F T P(k + 1) = QX(k) + F T K(k + 1)X(k + 1) (2.42) 

Da equação 2.42 obtém-se que: 

X(k+1)=FX(k)-Gz~ 1 G T K(k + 1 ) X(k + 1 ) (2.43) 

Da equação 2.43, anterior, isola-se a expressão para a matriz X(k+1). Assim, 
tem-se que: 


X(k+1) = [l + Gz~ 1 G T K(k +1 )]~ 1 FX(k) (2.44) 

onde lé a matriz identidade 

Levando-se a expressão de X(k+1) da equação 2.44 na equação 2.42 vem 
que: 


K(k)X(k) = QX(k)+ F T K(k + 1)[l + GZ 1 G T K(k+1)] 1 F X(k) 

Eliminando-se o vetor variáveis de estado X entre os membros da equação 
precedente vem: 

K(k) = Q+F T K(k + 1)[l + GZ~ 1 G T K(k + 1)]~ 1 F (2.45) 

A equação 2.45 é a equação de Riccati em sua forma discreta que realmente 
contém a relação iterativa possibilitando, através de uma implementação 
computacional, que se obtenha a matriz K. 

De posse da matriz K, através de um simples procedimento matemático, 
obtém-se a matriz P(k+1) que, levada na equação 2.39, fornecerá a matriz das 
variáveis de realimentação da equação U(k). 

Da equação 2.38 extrai-se a expressão da matriz P(k+1) da seguinte maneira: 

P(k + i) = [F T r 1 [P(k)-QX(k)] (2.46) 


Levando-se a expressão de P(k+1) da equação 2.46 na equação 2.39 vem 
que: 
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U(k) = -Z 1 G T [F T ] 1 [P(k)-QX(k)] 

Levando, agora, na equação anterior, a expressão de P(k) da equação 2.41, 
obtém se: 

U(k j^-Z^^ 7 [F 7 f 1 [K(k)X(k)-Q X(k jj^-Z - ^ 7 [F 7 f 1 [K(k)-Q] X(k) = r X(k) 

(2.47) 

A equação 2.47 reúne, a partir da determinação da matriz de realimentação T, 
os elementos que garantem o controle ótimo via realimentação de estado. 

Para voltar ao plano contínuo bastará utilizar a fórmula seguinte: 

r (t) = T (k) (2.48) 


2.6. Exemplo Prático 

O exemplo aqui apresentado se refere ao sistema de regulação de velocidade 
da UHE Emborcação (4 x 298 MW) localizada no Triângulo Mineiro. A figura 
2.1 mostra o diagrama em blocos no plano de Laplace do sistema de regulação 
de velocidade. 

Os resultados podem ser vistos nas tabelas de 2.1 a 2.10 que mostram uma 
situação em que o sistema de regulação de velocidade passou por uma 
instabilidade e depois com a estabilização. Isso pode ser visto e confirmado 
pela presença de autovalores com partes reais positivas. Ver tabela 2.8. Logo 
após, nas mesmas condições de instabilidade, foi aplicado o procedimento com 
a minimização de índices de desempenho com a aplicação da técnica de 
realimentação de estado, via equação de Riccati. 

Valores em Por Unidade da Figura 2.1 

Potência Base = 284 MW 
Frequência Elétrica Base = 377 rad/s 
bp = 0,05 pu/pu 
Td = 5,2 s 
TD = 52 s 

KLP = 4,1 sem dimensão 
TI = 0,685 s 
Ty = 0,2 s 
Ymáx = +1,00 pu 
Ymín = -1,00 pu 


35 



Tw = 1,55 s 
2H = 10,43 s 
D = 0,5 pu/pu 

T = 0,01 pu/rad em rede isolada ou off-line 
T =3,1 pu/rad em rede interligada ou on-line 
Cl : comutação em 60 s 
C2 : comutação em 50 s 
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Bloco de Realímentação de Estado 


Figura 2.1 - Sistema Regulador de Velocidade 
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Observar na tabela 2.10 de autovalores que o sistema, dentro das ponderações 
adotadas, tornou-se estável. As figuras 2.2 e 2.3 mostram as respostas do 
sistema para a s duas condições. No instante 50s houve a comutação entre as 
operações on-line para off-line na rede de transmissão. Para efeito de impacto, 
no instante 20 s houve uma variação de 0,05 pu de carga. Na figura 2.3, a 
entrada do controlador via realimentação de estado se deu no instante 110 s. 


Tabela 2.1 - Matriz "F" Discretizada 


0.9977 

0.0000 

0.0000 

0.0000 

0.0418 

0.0000 

-0.2702 

0.7295 

-0.0005 

0.0003 

0.5349 

0.0000 

-0.0070 

0.0418 

0.9512 

0.0000 

0.0139 

0.0000 

-0.0001 

0.0009 

0.0375 

0.9872 

0.0002 

0.0000 

0.0000 

0.0000 

-0.0019 

0.0010 

0.9995 

0.0000 

0.0000 

-0.0001 

-0.0035 

0.0018 

3.7691 

1.0000 


Tabela 2.2 - Matriz "G" Discretizada 

0.0418 

-0.3945 

- 0.0102 

- 0.0001 

0.0000 

0.0000 


Tabela 2.3 - Matriz "K" De Riccati 


6.7787 

-0.3736 

-2.4090 

-8.9119 

-108.0920 

0.1410 

-0.3736 

0.0489 

0.3090 

1.8252 

26.6118 

-0.0174 

-2.4090 

0.3090 

2.0467 

12.3740 

181.5715 

-0.1149 

-8.9118 

1.8252 

12.3740 

181.0163 

3188.3181 

1.8181 

-108.0903 

26.6115 

181.5700 

3188.3149 

57324.2734 

40.4868 

0.1410 

-0.0174 

-0.1149 

1.8181 

40.4862 

0.0984 


Tabela 2.4 - Matriz de Realimentação 


0.4502 -0.0370 -0.2378 -1.1787 -16.2569 0.0136 
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Tabela 2.5 - Matriz "Q" De Entrada 


0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0010 


Tabela 2.6 - Matriz "A" Continua Original 


-0.2320 

0.0000 

0.0000 

0.0000 

4.1860 

0.0000 

-31.5380 

-31.5380 

0.0000 

0.0000 

63.0770 

0.0000 

0.0000 

5.0000 

-5.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

3.8700 

-1.2900 

0.0000 

0.0000 

0.0000 

0.0000 

-0.1917 

0.0958 

-0.0479 

-0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

377.0000 

0.0000 


Tabela 2.7 - Matriz "B" Continua Original 

4.1860 

-45.3038 

0.0000 

0.0000 

0.0000 

0.0000 




Tabela 2.8 - 

Autovalores Do Sistema Original 

REAL= 

-31.6196 

IMAGINÁRIO 

= 0.0000 

REAL= 

-3.3840 

IMAGINARIO= 

0.8184 

REAL= 

-3.3840 

IMAGINARO 

-0.8184 

REAL= 

0.1910 

IMAGINARIO= 

0.6458 

REAL= 

0.1910 

IMAGINARIO= 

-0.6458 

REAL= 

-0.1024 

IMAGINARO 

0.0000 
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Tabela 2.9 - 

Nova Matriz 

"A " 

-2.1167 

0.1547 

0.9954 

4.9338 

72.2372 

-0.0571 

-11.1406 

-33.2127 

-10.7730 

-53.3974 

-673.4207 

0.6174 

0.0000 

5.0000 

-5.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

3.8700 

-1.2900 

0.0000 

0.0000 

0.0000 

0.0000 

-0.1917 

0.0958 

-0.0479 

-0.0010 

0.0000 

0.0000 

0.0000 

0.0000 

377.0000 

0.0000 


Tabela 2.10- Autovalores do Sistema com Realimentação 


REAL= 

-31.6454 

IMAGINARIO= 

0.0000 

REAL= 

-3.4230 

IMAGINARIO= 

0.7150 

REAL= 

-3.4230 

IMAGINARIO= 

-0.7150 

REAL= 

-1.2709 

IMAGINARIO= 

1.6084 

REAL= 

-1.2709 

IMAGINARIO= 

-1.6084 

REAL= 

-0.6340 

IMAGINARIO= 

0.0000 
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Figura 2.2 - Sistema sem Realimentação de Estado 


41 












0,2 
0,1 
0 

- 0,1 
- 0,2 

0 20 40 60 80 100 120 140 160 180 200 


= Curso do Servomotor Principal em pu 


Tempo em s 





0 20 40 60 80 100 120 140 160 180 200 


Figura 2.3 - Sistema com Realimentação de Estado 


42 



















O procedimento iterativo utilizado para a computação digital da matriz K da 
equação de Riccati teve seu sucesso garantido a partir do momento em que 
houve convergência dos valores sempre a partir de condições iniciais 
favoráveis. Algumas observações se fazem necessárias: 

a) Durante o processo de cálculo da matriz da planta no plano discreto, 
matriz F, deve-se assegurar que o passo de integração adotado seja no 
máximo 10% inferior a menor constante de tempo envolvida no processo. 

b) Cuidados especiais deverão ser tomados com subsistemas com laços de 
realimentação. 

c) A condição de estabilidade aqui adotada consiste em considerar, a partir 
do prévio conhecimento das equações fásicas na forma canônica de 
variáveis de estado, que um determinado sistema somente será declarado 
estável se, e somente se, todos os seus autovalores associados à matriz da 
planta do modelo contínuo não apresentarem partes reais positivas. 

d) A matriz de ponderação relativa às variáveis de entrada, matriz Z no plano 
discreto, se tomada como diagonal unitária será suficiente para o bom 
desenvolvimento do processo computacional. 

e) A matriz Q de ponderação relativa às variáveis de estado tem importância 
capital no desempenho do processo. Devem-se testar várias matrizes a 
começar por uma matriz diagonal, porém escolher e testar os valores antes. 

f) A matriz K resultada a partir da solução da equação de Riccati será 
sempre simétrica desde que a matriz de ponderação Q o seja. 

g) Partir sempre da condição em que a matriz Q de ponderação seja uma 
matriz definida positiva, no mínimo semi-definida positiva. 

h) Por fim saber que a precisão do modelo é imprescindível ao bom 
desempenho dessa proposta de controle. 
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2.7. Considerações Finais 


A técnica do controle ótimo por realimentação de estado apresenta algumas 
vantagens sobre outros tipos de sistemas quando se quer prover um sistema 
de controle da estabilidade operativa eficiente e de baixo custo. No caso 
apresentado houve uma ocorrência em que os disjuntores de 500 kV, referente 
à condição de operação on-line do sistema, abriram todos ficando a UHE 
Emborcação alimentando apenas o sistema regional de 138 kV o que 
ocasionou instabilidade e, por conseguinte, o desligamento da unidade gerando 
um tipo de shortage. Essa ocorrência foi retratada na simulação cujos 
resultados podem ser vistos nos gráficos das figuras 2.2 e 2.3. Percebe-se que 
em 1 min, com a entrada do controlador por realimentação de estado, o 
sistema tornou-se estável. A vantagem da instalação do controlador por 
realimentação de estado é que ele pode ser ajustado antecipadamente e 
mantido ligado ao sistema regulador de velocidade com a sua entrada em 
funcionamento controlada externamente por dispositivos tais como o de 
monitoramento do estado atual dos disjuntores de 500 kV, no caso. Há 
situações em que o monitoramento poderá ser desenvolvido em função de uma 
unidade de proteção direcional de supervisão da direção do fluxo de potência 
ativa elétrica em uma determinada linha de transmissão, por exemplo. 

A implementação de um dispositivo do tipo controlador por realimentação de 
estado, utilizando-se da tecnologia da eletrônica digital, é mais apropriada 
ainda em função de sua efetivação no campo uma vez que por ação da 
modelagem numérica é possível o acesso a todas as variáveis de estado 
envolvidas no processo. Outra facilidade a ser introduzida seria o processo de 
ajuste no campo o que exigiria apenas ação sobre potenciômetros quer sejam 
estes analógicos ou digitais. 

Em algumas instalações elétricas, quando há perda da condição de operação 
on-line, o que pode ser caracterizado pelo desligamento dos disjuntores de 500 
kV, por exemplo, já existem instalados sistemas suplementares conhecidos 
como SAGE, ou subsistema de alívio de geração, cuja função é reduzir, 
durante tais ocorrências, a geração das máquinas uma vez que não se dispõe 
de carga suficiente versus a geração despachada antes da ocorrência. Nesses 
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casos, como há uma redução brusca no valor do torque sincronizante, mesmo 
com ajustes para rede isolada, é possível que o regulador de velocidade não 
consiga manter as máquinas interligadas. Esse foi o caso mostrado na 
simulação apresentada nesse artigo. 
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CAPÍTULO 3 - FENÔMENO DE AUTOEXCITAÇÃO EM GERADORES DE 

ENERGIA ELÉTRICA 


RESUMO: O fenômeno de autoexcitação em máquinas síncronas tem por 
razão os efeitos de ressonância elétrica aos quais estes equipamentos ficam 
sujeitos em função da presença de cargas indutivas associadas a cargas 
capacitivas. O resultado deste fenômeno pode ser acompanhado de fortes 
sobretensões tanto no campo quanto no estator da máquina. Esse fenômeno 
também pode ser agravado durante as oscilações no sistema interligado uma 
vez que o processo de ressonância tem a ver com a frequência elétrica 
envolvida. Este artigo tem por objetivo apresentar uma análise sucinta do ponto 
de vista analítico para, posteriormente, calcular a faixa de risco operativo 
considerando o fenômeno da autoexcitação. Serão utilizadas as chamadas 
transformações ortogonais de Park, também conhecidas como transformações 
de Blondel [2], seu primeiro idealizador. Na realidade a transformação a ser 
utilizada será uma variante daquela que foi sugerida por Anderson [1], 


Palavras-chave: Máquina síncrona, regulação de tensão, sistema de 
excitação. 


3.1. Explanação Geral Acerca de Máquinas Síncronas 

Ao se iniciar o estudo relativo à máquina síncrona, pode-se apresentar de início 
uma comparação direta entre as máquinas síncronas e outras máquinas 
rotativas. 

Com relação às máquinas assíncronas ou máquinas de indução, outro tipo 
físico de máquina elétrica, pode-se afirmar que: 

a) as máquinas síncronas, quando operam como gerador, fornecem, ao mesmo 
tempo, potência ativa e potência reativa, já o gerador assíncrono fornece 
apenas potência ativa, continuando como no motor assíncrono, a consumir 
potência reativa; 

b) um gerador assíncrono consome em média, por característica, uma potência 
externa de aproximadamente 40% de sua potência nominal no processo de 
excitação o que inviabilizaria sua utilização já que o processo de excitação de 
um gerador síncrono exige apenas uma potência da ordem de 1% para menos 
de sua potência nominal. 

A máquina síncrona pode ser representada, esquematicamente, pela figura 3.1. 
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ENROLAMENTO 



Figura 3.1 - Esquema Básico de uma Máquina Síncrona com um Par de Pólos 

e Ligação do Estator em Estrela 


Pode-se notar pela figura 3.1 que existem dois eixos, um na direção do pólo e 
outro em quadratura com este. O eixo na direção do pólo é chamado de eixo 
direto ou apenas d, enquanto o outro é chamado de eixo em quadratura ou 
apenas q. Apenas lembrando que nesse caso, por questões didáticas, o estator 
se refere à parte fixa da máquina e o rotor à parte que se move com velocidade 
co em relação ao estator. Observar que este sistema de eixos d-q, idealizado, 
gira na mesma rotação do rotor, ou seja, com velocidade co. 

Como o torque mecânico é transferido através da rotação do rotor via eixo de 
acoplamento rotor/fonte de energia primária, quando em gerador síncrono, ou 
via eixo de acoplamento rotor/carga mecânica, quando em motor síncrono, 
existe um sincronismo entre o campo magnético de excitação e o torque 
mecânico que é desenvolvido neste eixo de acoplamento, ou seja, ambos 
deverão girar na mesma velocidade, portanto sincronizados em regime 
permanente. Esta é a razão da máquina se chamar máquina síncrona. 

Nota-se, também, que o sistema de eixos d-q gira sincronizado na freqüência 
de rotação do conjunto girante uma vez que se trata de eixos fictícios 
idealizados sobre o rotor da máquina. A idealização desse sistema de eixos 
vem ligada às facilidades de análise quando todas as grandezas, mesmo 
aquelas localizadas no estator, são referidas ao rotor, portanto aos eixos d-q. 
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Estas transformações são clássicas no estudo da máquina síncrona e são 
conhecidas como transformações de Biondel, primeiro estudioso que as 
idealizou. Logo após, Park, Doherty, Nickle, entre outros estudiosos, 
apresentaram trabalhos nesta área. Nesse livro essas transformações serão 
chamadas apenas de transformações d-q-0. 

O estudo da máquina síncrona, considerando-se as transformações de Biondel, 
será então desenvolvido de maneira fundamentada uma vez que sem esse 
conceito é difícil compreender perfeitamente o processo de regulação de 
tensão da máquina síncrona. 

3.1.1. Enrolamento de Campo 

O enrolamento de campo que envolve cada pólo cumpre o papel de, a partir de 
uma fonte de corrente contínua, que pode ser externa ou auto-alimentada, 
produzir o campo magnético indutor responsável pela excitação da máquina 
síncrona. 

Em uma máquina síncrona com um par de pólos, por exemplo, nota-se que o 
fluxo magnético gerado em um dos lados do eixo polar é norte enquanto que 
do lado oposto o fluxo magnético gerado é sul. Em outras palavras, trata-se de 
um grande eletroímã. O enrolamento em si é constituído por espiras na forma 
de barras de cobre eletrolítico numa disposição retangular que envolve cada 
pólo de modo concentrado. 

3.1.2. Estator 

O estator da máquina síncrona, na sua forma padrão, constitui a parte estática 
que é varrida pelo fluxo magnético contínuo gerado no rotor. É formado pela 
carcaça, pelo núcleo magnético e pelos enrolamentos de fase da máquina. 

O núcleo do estator, semelhantemente ao núcleo do rotor, é confeccionado de 
chapas de aço silicioso, ou seja, uma liga de aço com silício que apresenta alta 
permeabilidade magnética e baixas perdas, isoladas umas das outras e 
montadas segundo um alinhamento definido previamente através de pinos guia 
e presas por tirantes. O conjunto é chumbado na estrutura civil da instalação 
elétrica. 


49 



O enrolamento do estator é constituído de bobinas cujos condutores 
encontram-se alojados em ranhuras na sua parte interna. A estrutura final do 
anel magnético do estator termina em uma armação circular tendo na sua parte 
externa a carcaça e os trocadores de calor e em sua parte interna uma 
superfície dentada com as ranhuras que comportarão as bobinas. O estator da 
máquina síncrona também é conhecido como armadura da máquina. 

3.1.3. Sistema de Excitação 

Este sistema é responsável peia produção do fluxo magnético indutor a partir 
de uma fonte de corrente contínua podendo ser externa ou auto-alimentada 
que, no caso de grandes máquinas, se constitui de um controlador 
automatizado e que em pequenas máquinas pode ser composto de um sistema 
reostático controlado manualmente. Nas grandes instalações, esse sistema 
torna-se complexo e cumpre importantes ações que se refletem na própria 
estabilidade do sistema elétrico de potência, SEP. 

A tensão de corrente contínua de alimentação é levada ao enrolamento de 
campo através de escovas que se encontram em contato permanente com os 
anéis coletores. 

Além daquilo que foi descrito até aqui, é importante lembrar que as máquinas 
síncronas se diferem quanto à construção do rotor e, por conseguinte, do 
enrolamento de campo. 

3.2. Equações da Máquina Síncrona 

3.2.1. Geral 

A máquina síncrona caracteriza o elemento representativo do processo sobre o 
qual vão atuar os órgãos de comando e controle para a obtenção do produto 
final, qual seja, a potência elétrica despachada a partir da potência motriz de 
entrada fornecida pela turbina. 

Para o equacionamento da parte relativa à regulação de tensão, faz-se 
necessária a utilização da representação da máquina síncrona através de um 
equacionamento tal que venha a facilitar as análises e os estudos envolvidos. 
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O procedimento normal é a aplicação das chamadas transformações d-q-0 que 
busca a conversão de grandezas estatóricas em grandezas rotóricas através 
da transformação de grandezas estacionárias em grandezas pseudo- 
estacionárias referidas a um sistema de eixos fixado sobre o rotor e que gira na 
mesma velocidade deste, velocidade síncrona, fazendo com isto que estas 
grandezas independam da rotação. 

O fato consiste na criação de três eixos chamados d-q-0, ou mais 
propriamente, eixo direto, eixo em quadratura e eixo zero, sendo que os eixos d 
e q perfazem entre si um ângulo de 90 graus elétricos e se movimentando na 
velocidade síncrona, sendo o eixo zero um eixo estacionário simplesmente 
utilizado nas representações em que as condições desequilibradas são 
verificadas. 

O processo visa à projeção das grandezas dos eixos estatóricos, aqui 
designados por a-b-c , segundo uma matriz de transformação particular, sobre 
aqueles eixos d-q-0, do rotor, e a partir daí obter o equacionamento de toda a 
dinâmica da máquina síncrona. 

Com o estabelecimento do equacionamento básico relativo às indutâncias, 
segundo as influências próprias e mútuas dos enrolamentos envolvidos, serão 
equacionadas as funções de transferência das partes implicadas buscando-se 
com isto a definição do modelo da máquina síncrona. 

Em cima da modelagem e da análise, detalhadas previamente, será possível a 
realização de diagnósticos sobre o comportamento operativo da máquina 
síncrona sob determinadas condições. 

3.2.2. A Transformação d-q-0 

Como já foi citada anteriormente, a idéia da utilização da transformação de 
grandezas estacionárias em grandezas pseudo-estacionárias teve seu 
embasamento matemático formulado na utilização de uma matriz de 
transformação, um método específico de tratamento matemático para essa 
questão. 
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A transformação a ser aqui utilizada se trata da transformação específica 
proposta por W A Lewis, ver referência [8]. 

Primeiramente será representada uma máquina síncrona composta de um par 
de pólos conforme está mostrado na figura 3.1. 

Os eixos d e q giram com a velocidade síncrona do conjunto no sentido horário. 
O eixo d encontra-se em fase com a direção do pólo enquanto o eixo q 
encontra-se em fase, na direção, 90 graus elétricos atrasados do eixo d, 
considerando o sentido de giro. 

O diagrama fasorial para a máquina síncrona, anteriormente mostrado na figura 
3.1, é apresentado na figura 3.2 considerando-se apenas a fase a do estator. 

A equação básica da transformação d-q-0 é a seguinte: 


G dqo = P G abc 


(3.1) 


onde 


G: matriz das grandezas no sistema indicado; 


P: matriz da transformação. 


Sentido de Rotação 



Eixo d 



Figura 3.2 - Diagrama Vetorial da Máquina Síncrona para a Fase a 
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Na figura 3.2, considerar que o ângulo elétrico de giro do sistema seja dado 
pela seguinte expressão: 


G(t) = õ+u)t + ^ (3.2) 

onde co é a rotação ou frequência do conjunto girante, geralmente tomada 
como um valor fixo, síncrono e sob controle automático em função da presença 
do regulador de velocidade. 

Em relação à figura 3.2, deve-se notar que essa representação expressa a 
posição momentânea dos vetores de tensão da máquina síncrona, mormente o 
vetor de tensão terminal da fase a, vetor v a em relação ao vetor da tensão de 
excitação E a que se encontra sobre o eixo q. Deste modo, num dado momento 
em que a localização do rotor fica definida, o vetor da tensão terminal relativa à 
fase a está localizado exatamente a um ângulo 0 radianos elétricos atrás do 
eixo d cuja flecha positiva indica o sentido da fmm perpendicular ao plano de 
uma espira hipotética que concentraria toda a corrente gerada na fase a pelo 
processo de conversão eletrodinâmica da máquina elétrica. A equação 3.2 é 
bem clara nessa interpretação. 

Nota importante: A partir desse ponto, no sentido de dar mais clareza e 
simplicidade ao equacionamento que será um tanto quanto trabalhoso, nas 
designações de grandezas variantes no domínio temporal será suprimida a 
denotação (t) como indicação de dependência da variável tempo. 

Assim, a equação 002, de modo mais simplificado, tomará a seguinte forma: 


d = Õ+LOt + — 

2 

Considerando- se a matriz P, da transformação desejada, vem que: 



sinQ 


7 2 


72 


2 2 
cos (9 — n) cos (9 + —n) 
3 3 

2 2 
siri (9 — n) sin(9 + —n) 
3 3 


(3.3) 
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A transformação P, aqui utilizada, trata-se de uma transformação ortogonal, ou 
seja, aquela que obedece à seguinte operação quando a inversa da matriz 
original é exatamente igual a sua transposta: 

P~ 1 =P T (3.4) 

Portanto, para as correntes, para as tensões e para os fluxos magnéticos 
enlaçados, grandezas variantes no tempo e associadas à armadura da 
máquina síncrona e refletidas no sistema d-q-0, pode-se escrever, de forma 
matricial resumida, que: 


'dqo = P 'abc 

(3.5) 

v dqo = P v abc 

(3.6) 

'dqo = P A abc 

(3.7) 


A matriz inversa da matriz P de transformação é a seguinte: 


P 




cosQ sinQ 

2 2 

cos (6 —tt ) sin(6 — n) 
3 3 

2 2 

cos(9 + —n) sin(6 + —n) 
3 3 


(3.8) 


A transformação escolhida garante a invariância da potência em ambos os 
sistemas a-b-c e d-q-0 e para comprovar isso bastará considerar a expressão 
da potência instantânea, no domínio do tempo, deliberada por uma máquina 
síncrona operando como gerador. Assim, vem que, levando-se em conta as 
tensões e correntes instantâneas de cada uma das fases do estator da 
máquina síncrona, variantes no tempo: 

P= v a i a + v b i b + v c i c < 3 - 9 ) 


Onde: 

p: potência instantânea total gerada; 

v a ,b,c: tensões instantâneas de estator por fase; 

ia,b,c : correntes instantâneas de estator por fase. 
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De outra maneira, a equação 3.9 pode ser reescrita numa forma mais 
adequada para as operações matriciais futuras. Deste modo, tem-se que: 




(3.10) 


Notar que as tensões e as correntes instantâneas de estator da máquina 
síncrona se distribuem na forma de uma matriz coluna, assim como tais 
entidades são denominadas na teoria. Isso quer dizer que: 


e 


i v abci 




v abc 


(3.11) 



\.'abc\ 


'b 



'abc 


(3.12) 


Notar que na equação 010, a matriz coluna das tensões de fase teve de ser 
transposta para que a operação de multiplicação matricial tivesse consistência. 

A partir das equações 005 e 006, podem-se escrever as seguintes equações, 
respectivamente para as correntes e as tensões de fase da máquina síncrona: 


'abc p 'dqo e v abc p v dqo 

Levando-se as expressões anteriores, respectivamente para as correntes e 
para as tensões, na equação 3.10, vem que: 


P = 


P~ 1 v 


dqo 


77 


P~ 7 / 


dqo 


v 


dqo\ 


7-1 


dqo < a13 > 


Uma importante observação relativa à transformação via matriz P é que ela se 
trata de uma transformação ortogonal, valendo então a relação da equação 
3.13 anterior e as observações concernentes a esta. 
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A partir da matriz de transformação P, tem-se, para a equação 3.13 referente à 
potência total instantânea transferida pela máquina síncrona, que: 

P =v dqo. 'dqo = v dqo. 'dqo (3- 14 ) 

Através do cálculo matricial anterior fica comprovada a invariância da potência 
na transformação a-b-c - d-q-0. Essa propriedade será fundamental para o 
estudo da transformação G^^ = P G abc e vice versa. 

Uma observação interessante: quando se trabalha com sistemas trifásicos 
equilibrados, as projeções das grandezas sobre o eixo 0 tornam-se nulas uma 
vez que a existência desse eixo tem a ver apenas com a questão de um 
possível desequilíbrio entre as grandezas trifásicas. 

3.2.3. Equação Geral dos Fluxos Enlaçados 

Na figura 3.1, embora muito simples, podem-se verificar os enrolamentos que 
interferem no processo de trabalho de uma máquina síncrona. 

Tendo em vista se tratar de uma máquina síncrona trifásica, notam-se os 
enrolamentos do estator, em número de três, o enrolamento de campo e os 
enrolamentos amortecedores. 


Em todos estes enrolamentos aparecem os efeitos da indução própria e da 
indução mútua, uma vez que eles se interagem simultaneamente. 


De uma maneira geral, podem-se equacionar os efeitos dos fluxos enlaçados 
segundo a relação matricial seguinte. 


''a 


L aa 

L ab 

L ac 

Á b 


L ba 

L bb 

L bc 

*c 


L ca 

L cb 

L cc 



L Fa 

L Fb 

l Fc 

a d 


L Da 

L Db 

L Dc 

o 


L Qa 

L Qb 

L Qc 


aF L aD L aQ 'a 
L bF L bD L bQ 'b 
L cF J~cD L cQ 'c 

l ff l fd l fq 'F 
l df l dd l dq 'd 
l qf l qd l qq 'q 


onde 
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Aa,b,c,F,D,Q : fluxos instantâneos enlaçados nos enrolamentos indicados pelos 
índices 

L a ,b,c,F,D,Q : indutâncias de estator e rotor por fase 

ia,b,c,F,D,Q : correntes instantâneas de estator e rotor por fase 

Considerar para a designação da indutância Ly que: 

j=i : para a indutância própria do enrolamento considerado; 

j*i : para a indutância no enramamento i devido ao enrolamento j. 


3.2.4. Indutâncias Próprias e Mútuas do Estator 

Para obtenção das fórmulas das indutâncias próprias e mútuas do estator de 
uma máquina síncrona considerar os cálculos cujos desenvolvimentos podem 
ser vistos na referência [6]. Simplesmente, as expressões são as seguintes: 


L aa = L s + L m cos2e 


(3.16) 


onde 


2 ' s c/ + ' s g 


L s = N a 


2 J d~ J q 


L m N a 


Nas equações anteriores considerar: 

N a : número de espiras da fase a do estator; 


3d: permeância magnética do eixo d do rotor; 


3 q : permeância magnética do eixo q do rotor. 


Portanto Laa possui uma parcela síncrona Ls e uma parcela dependente da 
velocidade dupla L m cos29. Considerando as outras duas fases do estator da 

máquina síncrona, têm-se as seguintes expressões para as indutâncias 
próprias do estator: 
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( 3 . 17 ) 


L bb= L s + L m cos2 ( 9 -^-) 

L cc= L s + L m cos2 ( e+2 ^) 
Para a indutância mútua tem-se que: 

L ab = L ba = - M s- L m cos2 ( 6+ \ > 

onde 



(3.18) 


(3.19) 



Portanto Lab ou Lba possui uma parcela síncrona Ms e uma parcela dependente 
da velocidade dupla L m cos 2 (9 + ^). 

Utilizando da analogia com a equação 3.19 relativa à influência da fase a na 
fase b, considerando-se que as outras duas fases do estator da máquina 
síncrona devem receber tratamento matemático semelhante, já que o sistema é 
considerado trifásico e equilibrado, tem-se para estas fases as seguintes 
expressões para as indutâncias mútuas do estator: 

L bc= L cb = - M s- L m cos2 ( Q ( 3 ' 2 °) 

5tt 

L ca = L ac = - M s- L m c ° s2 < e + -^> < 3 ' 21 > 

3.2.5. Indutâncias Próprias e Mútuas do Rotor 

As indutâncias próprias e mútuas do rotor independerão do ângulo 0 de giro da 
máquina síncrona, referenciado na equação 002, uma vez que seus valores já 
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se encontram vinculados aos eixos d e q relativos ao sistema de eixos já em 
movimento na velocidade síncrona. 

Assim, para as indutâncias próprias do rotor, tem-se, considerando no rotor os 
enrolamentos amortecedores de eixo direto e de eixo em quadratura D e Q, 
respectivamente, e o enrolamento de campo apenas no eixo direto, 
enrolamento F, seguindo a designação mais comum na literatura. 

Para o enrolamento de campo tem-se que: 




Para o enrolamento amortecedor de eixo direto: 


L DD ~ L D 

Para o enrolamento amortecedor de eixo em quadratura: 

L QQ = L Q 


(3.22) 


(3.23) 


(3.24) 


Para as indutâncias mútuas tem-se que: 
para campo e o enrolamento amortecedor D: 

L FD = L DF = M R 

para campo e o enrolamento amortecedor Q: 

L FQ = l QF = 0 

para os enrolamentos amortecedores D e Q: 

L DQ = L QD = 0 


(3.25) 


(3.26) 


(3.27) 


Nota: nas equações 3.26 e 3.27, as indutâncias mútuas entre os enrolamentos 
D e Q e entre os enrolamentos F e Q, nesse sentido, são nulas, isso porque 
esses enrolamentos fictícios perfazem entre si um ângulo teórico de 90 Q . O 

r — > „ 

fluxo indutor da máquina síncrona é dado pela fórmula 0 = ] B-nds, uma 
grandeza escalar do cálculo vetorial onde o integrando é o produto escalar do 
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vetor indução magnética B peio vetor relativo à área ds=hds, um vetor 
normal à superfície S por onde o fluxo passa. Essa integral, para o presente 
caso, é nula uma vez que o ângulo envolvido vale 90 Q . 

3.2.6. Indutâncias Mútuas Entre Estator e Rotor 

Considerando que a fem gerada no enrolamento de campo apresenta apenas a 
componente de eixo direto, torna-se possível a obtenção das expressões para 
as indutâncias mútuas do estator com o enrolamento de campo. 

L aF = N f 2 ^d cosd (3.28) 

ou simplesmente 

L aF = L Fa = MpCOsQ (3.29) 

Analogamente, para as demais fases, tem-se que: 

L bF= L Fb = m F cos < 9 - > < 3 - 30 > 

L CF = L Fc = M F C0S < 6+ > < 3 - 31 > 

Para as indutâncias mútuas entre as fases do estator e o enrolamento 
amortecedor de eixo direto, tem-se, considerando a fase a e o enrolamento D: 

L a u = Ljj a = Mjj cosG (3.32) 

onde 


M D-N a 2 3 d 


Analogamente para as indutâncias mútuas entre as demais fases do estator e o 
enrolamento amortecedor de eixo direto, tem-se: 


L bD~ L Db 


Mjj cos(Q- 



(3.33) 
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( 3 . 34 ) 


L cD = L Dc = M D C0S ( d+ ^-) 

Para as indutâncias mútuas entre as fases do estator e o enrolamento 
amortecedor de eixo em quadratura, Q, tem-se um caso análogo ao anterior, 
porém com a projeção executada sobre o eixo q. 

Assim 


L aQ~ L Qa ~ M Q sin6 

(3.35) 

L bO= L Qb = M Q si ^- Z 3 > 

(3.36) 

L cQ = L Qc= M Q sin < B+ 2 3 > 

(3.37) 


onde 

M Q = w a% 

3.2.7. Equações dos Fluxos Enlaçados e Desenvolvimento 

Considerando-se a equação 3.15, relativa aos fluxos enlaçados envolvidos com 
a máquina síncrona, e em seus valores instantâneos, numa forma simbólica 
mais compacta, vem: 


^ a,b,c 

Á F,D,Q 

Multiplicando os membros da equação 3.38 pela matriz seguinte indicada na 
equação matricial obtém-se que: 

"P 0 
0 I 


A a,b,c 

K F,D,Q 


P 0 
0 I 


L EE l er 
l re l rr. 


I 


a, b,c 


F,D,Q 


(3.39) 



l ee 

l er 




, l re 

l rr_ 


'f,d,q 


(3.38) 


onde 

I é uma matriz unitária 3x3 e P a matriz da transformação d-q-0. 
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Promovendo-se as operações, inclusive a multiplicação do segundo membro 
da equação 040 pelo produto unitário matricial P _1 P=1, vem que: 


A o,d,q 

Á F,D,Q 


p 

0 ~ 

\ l ee 

l er 


p~ 1 

0 


~P 0 


^, 0,0 

0 

1 

l re 

l rr_ 


0 

1 


0 1 


Jf,d,q_ 


(3.40) 


E também: 


Á o,d,q 


~P 0 

l ee l er 


P ~ 1 0 


1 0 ,d,q 

Á FD,Q_ 


0 1 

l re l rr_ 


0 1 


_'f,d,q_ 


A equação matricial 3.41 encontra-se numa forma compactada. 
Descompactando-a obtém-se a seguinte equação: 


onde 


Á 0 


L o 

0 

0 

0 

0 

0 

Á d 


0 

L d 

0 

km f 

km d 

0 

Á q 


0 

0 

L q 

0 

0 

km q 

a f 


0 

km f 

0 

l f 

m r 

0 

a d 


0 

km d 

0 

M fí 

l d 

0 

1 

O 


0 

0 

km q 

0 

0 

l q 


'O 

'd 

'q 

'F 

'D 

'q 


(3.42) 



L 


d 


L s + M s 




L s + M s 



L 0 = L s~ 2M s 

3.2.8. Equações das Tensões 

O esquema do diagrama da figura 3.3 mostra os circuitos dos enrolamentos 
componentes de uma máquina síncrona real operando como gerador. 
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Considerando o esquema da figura 3.3, segue então a equação matricial 
compactada das tensões envolvidas na máquina síncrona, em seus valores 
instantâneos. O ponto sobre a grandeza indica derivada desta em relação ao 
tempo. 


v-v n + r i-A+O 


(3.43) 


De maneira expandida, a equação matricial anterior toma a seguinte forma: 


v. 


v 


b 

v c 

~ V F 

0 

0 


v 


r a 

0 

0 

0 

0 

0 

0 

r b 

0 

0 

0 

0 

0 

0 

r c 

0 

0 

0 

0 

0 

0 

r F 

0 

0 

0 

0 

0 

0 

r D 

0 

0 

0 

0 

0 

0 

r c 



'a 


'b 


i c 


[F 


'D 


_'q_ 


a 

; 


c 

V 

'o 


= o 


(3.44) 




V" 

u L 


bb 


n 


3 Ln '^Lcc 


u,- aa 

L bb 


■ a 




cc 


n 


1 n 


•b 


■c 


a 


V C 


n 



Figura 3.3 - Enrolamentos Reais da Máquina Síncrona 
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Na figura 3.3 considerar, além das variáveis e parâmetros conhecidos e já 
definidos anteriormente, a seguinte legenda: 

r a a,bb,cc,F,D,Q : resistências dos enrolamentos; 

vf : tensão instantânea de campo; 

i n : corrente instantânea de neutro 

r n : resistência de neutro do estator; 

L n : indutância de neutro do estator. 

Para determinação da matriz v n relativa às correntes envolvidas tem-se: 


1 

cr 

cr 

cr 

i_ 


ia 


i-n Ln Ln 

c 

c 

c: 


ib 

— 

Ln Ln Ln 

i 

c: 

cr 

cr 

_i 


Jc_ 


Ln Ln Ln_ 


ia 

ib 

ic 


(3.45) 


De maneira mais resumida, a equação 3.45 toma a seguinte forma mais 
compacta: 


V„ = -RJ 


n 


n 


'a, b, c L n ' a > b’ c 


(3.46) 


De maneira compacta, a equação 3.43 poderá ser reescrita da seguinte forma: 


v a,b,c 
V F, D,Q 


v 


R 


a,b,c 

0 


0 

R F,D,Q 


’a,b,c 

'f, d,q 


ha,b,c 

hF,D,Q 


= 0 


(3.47) 


Aplicando-se à equação 3.47 a transformação d-q-0 pela matriz P, sobre as 
variáveis das vem que: 


~P 0 


v a,b,c 


v 0,d,q 

0 1 


V F, D,Q 


V F, D,Q 


(3.48) 


Considerando, a partir daqui, toda a dedução algébrica que pode ser vista nas 
referências [1] e [6], obtém-se a seguinte equação matricial das tensões 
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instantâneas em função das correntes instantâneas relativas aos enrolamentos 
da máquina síncrona. 


^0 


r a +3r n 

0 

0 

0 

0 

0 

'o 

v d 


0 

r b 

wL q 

0 

0 

lo KM q 

'd 

V Q 
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0 

~ WL d 

r c 

- CúMfr 

— üjKM d 

0 

'q 

- V F 

0 

0 

0 

r F 

0 

0 

'f 

V D = ° 


0 

0 

0 

0 

r D 

0 

'd 

V Q = 0 


0 

0 

0 

0 

0 

r Q 

'q 


L o +3L n 

0 

0 

0 

0 

0 

0 

L d 

0 

km f 

km d 

0 

0 

0 

L q 

0 

0 

km q 

0 

km f 

0 

l f 

M Ft 

0 

0 

km d 

0 

m r 

l d 

0 

0 

0 

km q 

0 

0 

l q 


i o 
id 







iF 
i D 
iQ 


(3.49) 


Utilizando a equação 3.49, pode-se representar o circuito equivalente para uma 
máquina síncrona, como se vê na figura 3.4, com r a = rb = r c = r , uma vez que 
todo projetista de máquina síncrona persegue esta condição. 

Pela figura 3.4, pode-se perceber, no circuito dos eixos d e q, a presença de 
duas tensões que são proporcionais ao movimento atual de rotação do 
conjunto e que são chamadas tensões de rotação. 

Portanto, as relações obtidas não são lineares, mas, por aproximação, 
tomando-se a velocidade de rotação co como constante, uma vez que a 
máquina síncrona dispõe um regulador de velocidade, o problema fica 
resolvido obtendo-se a linearização procurada. Lembrar que é mais prático e 
conveniente se trabalhar com equações lineares. 
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Figura 3.4 - Circuito Equivalente da Máquina Síncrona 

3.2.9. Valores em Por Unidade no Estudo da Máquina Síncrona 

Cabe aqui uma observação, para a utilização das equações da máquina 
síncrona em por unidade há necessidade de se proceder à normalização das 
grandezas envolvidas no processo. 

Para a normalização das grandezas utilizadas deverá ser adotado um 
referencial baseado no rotor, proposição feita por alguns autores apesar da 
escolha depender de certos critérios. 

Os fluxos enlaçados mútuos dependem da parcela mútua das indutâncias que 
compõem o sistema magnético da máquina síncrona. 

Considerando as parcelas do fluxo mútuo e do fluxo de dispersão, por 
definição, pode-se escrever a seguinte expressão para o valor total da 
indutância de eixo direto da máquina síncrona, primeiro caso a ser visto e a ser 
seguido daqui para frente neste tópico. 


L d ~ L md + L dd 


(3.50) 
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onde 


Lmd : indutância mútua útil referente ao eixo d; 


Lmd: indutância de dispersão útil referente ao eixo d. 

A indutância de dispersão tem a ver com o fluxo perdido no entreferro. 

Assim, ficam definidas, para todas as indutâncias em jogo, as relações 
seguintes referentes às parcelas úteis de cada uma em seu respectivo 
enrolamento conforme é mostrado na figura 3.4: 


md = L d - L dd 
Lmq = ~ L dq 

L mF = L F~ L dF 
L mD = L D~ L dD 
L mQ = L Q~ L dQ 


(3.51) 


As relações elementares para um determinado valor de fluxo enlaçado Ab, aqui 
tomado como um valor base que construirá a base de valores a ser escolhida 
no processo de normalização envolvendo outras grandezas tais como a tensão 
base Vb, a indutância base Lb , a corrente base Ib e rotação base cob, são as 
seguintes: 


Á B ~ L B 'b 


(3.52) 



V, 


B 


CO 


B 


~ l b'b 


(3.53) 


com 


uj b - 


l B 


(3.54) 


Notar que a equação 3.54 define uma grandeza que não deve ser confundida, 
trata-se do tempo base que é exatamente o inverso da velocidade angular 
elétrica, o que para o sistema elétrico nacional, SIN, tem seu valor aproximado 
para u)q = 2tt60 = 377rad/s 


67 



As seguintes premissas serão obedecidas: 


S B =V B J B potência para sistemas monofásicos; 

S B = J3V B I B -> potência total para sistemas trifásicos. 
onde 

Sb : MVA monofásico de base para sistemas monofásicos; 

Sb: MVA trifásico de base para sistemas trifásicos. 

Podem-se então definir, a partir de certas condições, as grandezas básicas do 
rotor. 

Assim, tem-se que: 


Á md = L md'B = KM F l FB = KM D 'üB 

Á mF = KM F 'b = L mF'FB = M R 'üB 
Á mD = KM D I b = M r I fb = L mD l DB 

Á mq = L mq'B = KM q'qB 

Á mQ = KM Q 'b = L mQ'QB 


(3.55) 


Nas equações 3.55, a ideia foi obter uma escolha de correntes base que, a 
partir da corrente base circulante por um determinado enrolamento, através de 
sua indutância própria útil, obter uma corrente induzida base nos outros 
enrolamentos por esse primeiro enrolamento influenciado por intermédio das 
indutâncias mútuas. 


EXEMPLO 3.1 


Utilizando-se das equações 3.55, provar que as indutâncias mútuas em pu 
KMFue KMdu são iguais. 

Solução: 

Das equações 3.55, podem-se escrever as seguintes relações: 
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Á md ~ L md l B ~ KM F 1 FB ~ KM D 1 DB 


Tomando-se a última igualdade tem-se: 


KM F 1 FB ~ KM D ! DB 


km F _ 1 DB 
KM D 1 FB 


(3.56) 


Porém sabe-se ainda da equação 3.53 que: 


km FB ~ 


V, 


B 


U B 1 FB 


'fb ~ 


V t 


B 


U B KM FB 


(3.57) 


KM 


V B 

DB U B 1 DB 


I 


V B 

DB u) B KM db 


(3.58) 


Levando na equação 3.56 as expressões encontradas nas equações 3.57 e 
3.58 para as correntes de base Ifb e Idb, obtém-se: 




B 


KM 


F u b km db km fb km f KM 


KM 


D 


V B 

üj b km fb 


KM 


DB 


KM 


FB 


KM FII =KM nil 

KM DB FU DU 


Por analogia ao resultado obtido anteriormente, pode-se generalizar para os 
circuitos dos eixos d e q da máquina síncrona que: 


L mdu ~ L mFu ~ L mDu ~ KM FU ~ KM DU ~ M RU 


(3.59) 


L mqu ~ L mQu ~ KM QU 

As equações 3.59 são de extrema importância quando da elaboração dos 
circuitos equivalentes de eixos d e q, em pu, da máquina síncrona. 

Nota a respeito do exemplo 1: 

Observar sempre que, na definição dos valores base para as indutâncias KMfb 
e KMdb, nas equações 3.57 e 3.58, a tensão base se refere ao enrolamento em 
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questão, porém a corrente base se refere ao enrolamento que criou o fluxo 
magnético. Essa proposta de base de valores é assim adotada para que os 
valores em pu das indutâncias mútuas sejam iguais, uma particularidade 
opcional visando à obtenção de facilidades para futuros cálculos. 


EXEMPLO 3.2 


Comprovar as equações 3.55 a partir de uma condição genérica. 

Solução: 

Para a comprovação propriamente dita, será utilizado nas variáveis um 
subscrito, apenas nesse exemplo, para que fique bem indicado a qual dos 
circuitos equivalentes d, D, F, q, Q o cálculo se refere. 

Imaginando uma escolha genérica para as correntes de base circulantes, tem- 
se o seguinte raciocínio. 

a) enrolamento d 

Os fluxos enlaçados vistos do enrolamento d serão os seguintes então: 


^ md d = L md l B 
Á mF d = L mF'FB 
Á mD d = L mD'DB 


(3.60) 


Para os três fluxos anteriores, a tensão base será sempre a tensão referente 
ao eixo d, ou seja, Vb, já que o efeito se dá no enrolamento de estator. 
Tomando-se os valores em pu para cada um dos três fluxos enlaçados da 
equação 3.60 vem que: 


^ Á md d = L md [b = L md 

mddu Á md dB L mdB'B L mdB 

A indutância base referente ao eixo d é dada por: 


(3.61) 
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(3.62) 



b) enrolamento q 

Estendendo o raciocínio para o enrolamento q vem: 

'b L mq = 'qB KM Q = A mq ( 3 - 68 ) 

c) enrolamento F 

Para o enrolamento de campo, em termos de fluxo magnético enlaçado, 
considerando-se a circulação da corrente base desse enrolamento, tem-se que: 

Á mF F = L mF l FB 
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Tomando o valor em pu para fluxo magnético enlaçado desse enrolamento, 
obtém-se: 


i _ L mF'FB 

mF Fu I l 

Fu L mFB'FB 


(3.69) 


A indutância base referente ao eixo F é dada por: 


L 


= V FB 


(3.70) 


Com a expressão anterior, a equação 3.69 pode ser assim reescrita: 


_ u B l FB 

ml ~Fu V FB mF 


(3.71) 


Por analogia, os fluxos enlaçados para AmdFu e AmDdu serão: 


(Ú Q /, 


777 d, 


B'B 


Fu V, 


KMr 


FB 


(3.72) 


mDi 


u) B l D 


Fu V t 


M 


FB 


R 


(3.73) 


Também aqui, para que, por suposição, os fluxos enlaçados em pu, relativos 
ao eixo em questão, sejam iguais, as seguintes relações deverão ser 
cumpridas: 

Á md Fu = A mF Fu = A mD Fu ( 3 - 74 ) 

Deste modo também, tem-se então que: 

'FB L mF = 1 B KM F = 1 DB M R = Á mF ( 375 ) 


d) enrolamentos D e Q 

Por analogia, podem-se escrever as seguintes equações para os fluxos 
enlaçados nesses enrolamentos: 
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'dB L mD ~ 'b KM D ~ 'fB M R~ A mD 


( 3 . 76 ) 


'qB L mQ ~ 'b KM Q ~ Á mQ 


(3.77) 


Ficando, deste modo, comprovadas as equações 3.65. 

Nota a respeito do exemplo 2: 

Multiplicando-se as equações 3.67, 3.75 e 3.76, relativas ao eixo d, 
respectivamente, por Ib, Ifb e Idb, vem que: 


Á md l B ~ L md l B ~ KM F 'FB 'b ~ KM D 1 DB 'b 


(3.78) 


A m F 'FB ~ KM F 'b 'FB ~ L mF 'FB ~ M R 1 FB 1 DB 


(3.79) 


e 


Á mD 1 DB ~ KM D 1 B ! DB ~ M R *FB 1 DB ~ L mD^DB ( 3 ■ 80 ) 

Do mesmo modo para o eixo q, multiplicando-se as equações 3.68 e 3.77, 
respectivamente, por Ib e Iqb, vem que: 

W B= L mq'B 2 = KM q'qb'b ( 381 ) 


e 


A mQ'QB ~ KM Q 'b 'QB ~ L mQ { QB ( 3 - 82 ) 

Por uma imposição do problema, as equações 078, 079 e 080, relativas ao eixo 
d, são iguais, logo vale a seguinte relação primordial: 

L md'B 2 = L mF'FB 2 = L mD , DB 2 = KM F 1 FB 'b = KM D 'üB 'b = M R 'FB 1 DB 

(3.83) 
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Do mesmo modo para o eixo q, as equações 3.81 e 3.82 são iguais, logo vaie a 
outra relação primordial: 


-mq'B 2 ~ L mQ'QB 2 ~ KM q'qB 'b 


(3.84) 


As equações 3.78, 3.79, 3.80, 3.81 e 3.82, chamadas aqui de equações 
primordiais, são relações que confirmam que as potências desenvolvidas nos 
diferentes enrolamentos dos eixos d e q da máquina síncrona, em termos dos 
valores base escolhidos nesse contexto, são iguais. 

Na realidade, as equações 3.55 são uma consequência da condição imposta 
relativa à conservação da energia entre os enrolamentos dos eixos d e q da 
máquina síncrona. 

Considerando as equações 3.83 e 3.84, além de se manter a mesma potência 
base de valores, podem-se escrever, para os três enrolamentos do rotor, as 
seguintes expressões: 

! 1 FB_ = Xb_ = \ L md _ L md _ KM F _ KM P _ K (3 Q5) 
'b V FB v L mF KM F L mF M fí F 


'_DB = ^B_ = \ L md _ L md _ KM P _ KM F _ K (3 86) 
'b V DB \ L m D KM D L mD M R° 


'QB 

'b 


V, 


B 


V, 


QB V 


L 


mq L mq 


KM, 


Q 


L 


mQ 


KM, 


L 


= K, 


Q mQ 


Q 


(3.87) 


Assim, as variáveis de corrente de base e de tensão de base são escolhidas de 
modo que em qualquer enrolamento a potência base seja sempre a mesma. 
Portanto, conhecendo-se Lb e Rb, das bases do estator, vem: 


r FB _ 1 R DB _ 1 

r b k f 2 ’ r b k d 2 


R, 


QB 


R 


B 


1 



(3.88) 
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1 


1 


(3.89) 



Kt 



K 


D 


L QB _ 1 

l b k q 2 


E para as indutâncias mútuas será utilizada uma propriedade extraída da 
equação 3.90, senão veja-se. 


L md 

km f 


KMl 


L 


mF 



^md^~mF 


(3.90) 


Nota: nessa propriedade, nota-se que, considerando o efeito mútuo entre dois 
circuitos magneticamente conectados, o produto das indutâncias próprias úteis 
de cada circuito elétrico é igual ao quadrado da indutância mútua verificada 
entre ambos. 


Deste modo, por conseqüência, levando-se em conta uma base de valores 
previamente definida, o quadrado do valor da indutância mútua base entre dois 
circuitos magneticamente conectados será igual ao produto dos valores das 
indutâncias próprias base de cada um dos circuitos envolvidos. 

Assim: 


M 


FB ~ L B L F B 


M FB l fb 


L 


B 


M 


(3.91) 


FB 


L 


Mas, da equação 085, tem-se que mQÍ = K F em valores atuais o que permite 


KM 


escrever em valores base que: 


Mas, da equação 3.86, tem-se que 
permite escrever em valores base que: 


M 


FB 


-B 


-md 


KM 


= K 


D 


D 



(3.92) 


em valores atuais o que 


m DB _ J_ 
L B k d 


(3.93) 
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= Kq em valores atuais o que permite 


Mas, da equação 087, tem-se que-— 

KM q 

escrever em valores base que: 


m QB _ 1 
L B k q 


(3.94) 


Da equação 3.85 tem-se que: 


km d 

m r 



L/nd 

L md 


km d 

m r 



m r 

L md 


1 KM 


D 


K F L md 


(3.95) 


Da equação 3.86 tem-se que: 


km d 

L md 



(3.96) 


Logo, a partir da equação 3.95, obtém-se para a equação 3.96, que: 


M R 

L md 


K F K D 


(3.97) 


Em termos de relação base, a partir da equação 3.97, vem que: 


M RB _ 1 

L B k f k d 


(3.98) 


A equação 098 completa o conjunto de equações envolvendo os valores base 
fornecendo o caminho para o cálculo das demais relações. 


3.2.10. Equações em Por Unidade 

Colocando a equação matricial 3.49, que se encontra em valores atuais, em 
por unidade obtém-se que: 
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1 QB lQu 

(3.99) 


Analisando-se cada uma das tensões do vetor tensão, em separado, iniciando- 
se peio eixo d vem que: 


v du + r ^ i du + “ L q^'qu 


IqB . Ir * 

+ ujKMq —— i Qu + L rí — i d u+ 


Vi 


B 


Vi 


B 


+ KM o — i Fu + KM d '-PP- iDu = 0 


(3.100) 


Vi 


B 


Vi 


B 


Sabe-se que 


1 

l B ^B 


Ir 

—, assim a equação 3.100 se transforma na equação 


3.101. 


v 


+ 


du + r u'du 


Ldu 'idu + 


CO 


B 


+ U) u L qu i qu + Lü u KM Qu i Qu 
KMr 


KM Fu + 
i Fu+ 


'Du 


ÜJ 


B 


(JÜ 


I Du = o 


B 


(3.101) 


Analogamente, para o eixo q, vem: 
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(3.102) 



r Qu'Qu 


KM Qu * L Qu * 

+- 'qu+ - iQu = O 

ÜJ g 6Ü g 


(3.106) 


Colocando as equações 3.101, 3.102, 3.103, 3.104, 3.105 e 3.106 na forma 
matricial completa, e normalizando o tempo nas parcelas que apresentam 
derivadas, obtém-se que: 
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Figura 3.5 - Circuito PU Equivalente de Correntes da Máquina Síncrona para 

(a) Eixo d, (b) Eixo q 
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= 0 


(3.107) 


3.2.11. Circuito Equivalente da Máquina Síncrona 

Como condição a ser obedecida nos desenvolvimentos que se seguem, deve- 
se levar em conta que o sistema é equilibrado eletricamente com io = 0. A base 
de valores escolhida é aquela descrita no item 3.2.10, precedente, cuja escolha 
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partiu de certas premissas a fim de se fazer com que as indutâncias mútuas 
fossem iguais em seus valores em pu. 

Os fluxos enlaçados da máquina síncrona, assim como se viu no item 3.2.9, 
podem ser relacionados, em pu, com os fluxos mútuos e os fluxos de 
dispersão, respectivamente, de modo que se pode escrever de maneira 
genérica que: 


A Tu- A du + A mu 


(3.108) 


ou seja, o fluxo enlaçado total At , de um determinado circuito magnético, é a 
soma dos fluxos de dispersão Ad e mútuo Am, desse mesmo circuito. 

Assim, tem-se que: 


A Tdu = A ddu + A mdu = L ddu'du + L mdu i du + KM Fu'Fu + KM D u'D u 
A Tqu = A dqu + A mqu = L dqu'qu + L mqJqu + KM Qu'Qu 
A TFu = A dFu + A mFu = L dFu'Fu + L mFu'Fu + KM Fu'du + M Ru'Du (3 ' 109) 
A TDu = A dDu + Á mCH = L dDu'Du + L m D dD u + M Ru'Fu + KM Du'd u 
A TQu = A dQu + A mQu = L dQu'Qu + L mQu i Qu + KM Qu'qu 


Das equações 3.59, sabe-se que: 

L mdu = L mFu = L mDu = KM FU = KM DU = M RU 
L mqu =L mQu = KM QU 

As duas equações anteriores permitem que se escrevam as seguintes relações 
em pu: 


L mdu ~ KM Fu - KM Du ~ M Ru - L id u 


(3.110) 


L mqu ~ L mQu ~ KM Qu ~ L iqu 


(3.111) 


E criando os fluxos enlaçados associados a cada uma das indutâncias gerais 
das equações 3.110 e 3.111 vem que: 


A i du ~ L i du ( 'du + 'Fu + 'Du> 


(3.112) 
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A i qu L i qu ('qu + 'Qu ) 


( 3 . 113 ) 


O significado físico das equações 3.112 e 3.113 pode ser representado 
considerando-se cada um dos eixos d e q como está mostrado na figura 3.5. 

Para o circuito equivalente completo, devem-se levar em conta as equações 
em tensão 3.107, já considerada a condição de se estar trabalhando com um 
sistema equilibrado trifásico com as resistências das fases do estator tomadas 
com o mesmo valor, assim como se tem procedido até aqui. 

Para o eixo d pode-se escrever que: 

v du + r u i du + -^ B A du+0} u A qu = 0 (3.114) 

Da primeira equação das equações 3.109 extrai-se a expressão para Adu total, 
o que leva a seguinte expressão: 


Á Tdu ~ L ddu 'du + L mdu 'du + KM Fu 'Fu + KM Du 'Du (3 - 115 ) 
Assim, a expressão da tensão 3.114 para o eixo d tomará a seguinte forma: 

v du + r u 'du + ~~ L ddu ' du+ ÜJ u Á qu + ('du + 'Fu + 'W = 0(3 ' 116 ) 

u) B v üj b 

De maneira similar, chega-se às equações de tensão para os demais 
enrolamentos. 

Para o enrolamento de eixo q tem-se que: 


v qu + r u 'qu + — L dqu 'qu~ UJ u A du + 
Para o enrolamento de campo tem-se que: 


L iqu 

üü b 


( i qu +i Qu)=° (3- 117 ) 


- v 


Fu + r Fu 'Fu 


u 


L dFu' Fu+ 

B 


('du +l Fu + 'Du)-° (3-118) 
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Para os enrolamentos amortecedores D e Q têm-se: 


v Du ~ 0 Du 'Du + 


— L dDu ’Du + (l du + 'Fu +1 Du) ~° í 3 - 119 ) 

U B U B 


v Qu - 0 r Qu j Qu + 


L dQu 'Qu 


+ -^< l qu +i Qu>-° 


(3.120) 


As figuras 3.6 e 3.7 mostram os circuitos equivalentes de eixo d e eixo q, 
completos em tensão, respectivamente. 



Figura 3.6 - Circuito PU Equivalente de Tensão da Máquina Síncrona para o 

Eixo d 


r Qu L d Qu r u 



Q) À d u 


Figura 3.7 - Circuito PU Equivalente de Tensão da Máquina Síncrona para o 

Eixo q 
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3.2.12. Autoexcitação em Máquinas Síncronas 

O fenômeno da autoexcitação em máquinas síncronas é um processo que 
pode ocorrer quando a máquina estiver conectada a um circuito elétrico de 
característica capacitiva. Deste modo, qualquer tipo de dispositivo capacitivo 
que esteja acoplado à máquina síncrona pode vir a provocar esse fenômeno. 

Para análise, será suposto o sistema mostrado na figura 3.8. 

Máquina C arga 



Figura 3.8 - Circuito Máquina Síncrona Conectada a uma Carga Capacitiva 



Figura 3.9 - Diagrama Vetorial para o Circuito Máquina Síncrona Conectada a 

uma Carga Capacitiva 
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No diagrama da figura 3.9 podem-se observar as grandezas envolvidas no 
processo sendo ii a corrente circulante, quando o circuito aposentava apenas 
uma determinada carga indutiva, e sendo i c a corrente capacitiva adiantada da 
tensão terminal v após a entrada da carga capacitiva. 

O seguinte equacionamento básico é possível. 

v = v cos(-6) = v cos 6 

d 

v q = v sin(-6 ) = -v sin 6 


Com 


/ =/' sin(n/2-Q 0 )= i cosQ 0 cosQ 0 = i li 

d c c d I c 

i = i r cos(n/2-6 0 )= i sin6 0 sind 0 = i li 

q u i q / c 

Considerando-se apenas a presença da carga capacitiva vem que: 


v 


v = ~j X c i c — = ~j X t 


Assim, as tensões Vd e v q tomam a seguinte forma: 


/ 

v =-' / < JL > = -(-)! =-(~iX c )i =X c i 

d I I q q d 

c c 


I 


V , 


^ =v(^) = (-)(i) = ~jX c (i)=X c i 

q I I d d q 

c c 


(3.121) 


onde |X c |=1/coCe y' = V^7. 


Considerando o fluxo enlaçado pelo eixo d, através das equações 3.109, 
desconsiderando o enrolamento amortecedor, tem-se que: 


Á du - L du'du + KM Fu'Fu 


(3.122) 


Para o fluxo enlaçado pelo campo, através das equações 3.109, 
desconsiderando o enrolamento amortecedor, tem-se que: 


Á Fu - L Fu'Fu + KM Fu'du 


(3.123) 
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A partir das relações primordiais anteriores, a equação 3.124 tomará o seguinte 
aspecto: 

e qu = KM Fu'Fu + L du 'du ~ L 'du 'du (3 ■ 125) 

Considerando-se agora a equação 3.123, a equação 3.125 poderá ser reescrita 
da seguinte forma: 

e qu = A du~ L du 'du 

Passando a equação anterior para valores atuais vem que: 

e q = U)A d~ U)L 'd i d (3 - 126 ) 

Das equações 3.101 e 3.102 obtêm-se as seguintes expressões para as 
tensões Vd e v q , desprezando-se as resistências envolvidas. 
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(3.127) 


V = ~CUA n 
d q 

v =(joA d 
q 


Das equações 3.121, a partir das equações 3.127, conclui-se que: 

OJ Á d — X c Í d 
CU Ag = — X c íq 

Assim, a equação 3.126 tomará a seguinte forma: 

e q = X c i d -wU d id 

O que, finalmente, conduz ao seguinte valor da corrente id: 

. _ cu C 

' d ~ 1-cu 2 L' d C Gq 


(3.128) 


(3.129) 


(3.130) 


Se o mesmo raciocínio for levado em consideração para um determinado 
enrolamento de eixo q, onde todas as condições forem respeitadas, a seguinte 
expressão para a corrente i q será obtida: 


'q 


cu C 


cu 2 L'C 


(3.131) 


A equação 3.104 fornece a expressão em por unidade da tensão de campo, 
onde a resistência de campo é considerada, o que permite que se escreva a 
seguinte expressão, considerando-se o fluxo através do campo da máquina 
síncrona: 


v Fu~ r Fu 1 Fu~ - 0 (3.132) 

Na equação 3.132, mais precisamente no termo da derivada do fluxo, a parcela 
1 /coB foi incorporada à variável \ Fu , o que garante estar a variável tempo sendo 
dada em segundos. 

O valor do fluxo enlaçado no campo da máquina síncrona será dado por: 
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Á Fu ~ r Fu ío( /- ~ 'Fu )dt 
Fu 


(3.133) 


KM Fu 

Multiplicando a equação 133, membro a membro, pela relação — vem 

L Fu 

que: 


KMc, i r r— f Ur,, 

rU a _ hU ff/í^i /7 hu 


I Fu I 
l Fu l Fu 


lo < KM Fu r - km Fu 'Fu > dt 
Fu 


Considerando as relações primordiais e sabendo-se da referência [1] que 

Lp, i 

T' do = - L ± L é a constante de tempo transitória de eixo direto a circuito aberto 
r 

Fu 

do estator da máquina síncrona , vem que: 


ío (e fd u ~ KM Fu 'Fu ) dt 

' do 


'qu 


(3.134) 


KMp u ifz u , da equação 

Trazendo na equação anterior 3.134 a expressão para 
3.125, obtém-se a seguinte expressão mais geral: 

V = T-íoíOfdu-lOqu-tLxlu-LduWW < 3 - 135 > 

1 do 

Passando a equação 135 para valores atuais obtém-se que: 

e q = -F\' 0 {e fd -[e -u>(L d -U d )i d l}dt 

1 do 

Ou ainda: 

0 . Q Q 

- í = yi{—-l- í -(L d - L ' d )i d ])dt (3.136) 

LO T d0 J(J (O (O U u u 

Trazendo na equação 136 a expressão da corrente id, da equação 3.130, tem- 
se que: 
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A equação 3.137 traz, finalmente, a função de transferência entre as tensões 
de campo etd e a tensão de excitação e q quando do carregamento da máquina 
síncrona por cargas capacitivas, considerando-se o eixo q. Trata-se de uma 
função de transferência de primeira ordem com um pólo real cujo valor 
depende dos parâmetros envolvidos no processo. 

É evidente que a combinação de valores entre os parâmetros pode designar 
um pólo positivo ou um pólo negativo o que terá uma influência importante na 
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questão da presença ou não de sobretensão no campo. A questão torna-se 
crítica quando o valor do pólo for negativo o que denota pela equação 137 que 
a resposta no domínio do tempo do sistema será regida por uma exponencial 
crescente o que representará, fisicamente, a presença de sobretensões 
elevadíssimas. 

Para o caso do eixo d, por analogia aos cálculos efetuados até aqui, a seguinte 
função de transferência poderá ser obtida: 


'd(S) 


T' qo S 


1 -co z L n C 

S + - f - 

T qo (1-(*> L’ C) 


'WS) 


(3.138) 


3.2.13. Análise da Autoexcitação em Máquinas Síncronas 


Pelo que foi determinado no item anterior é possível estabelecer as faixas para 
os limites de valores dos parâmetros da máquina síncrona em função da 
questão da sobretensão originária no fenômeno da autoexcitação provocado 
pela interação entre os parâmetros de máquina versus os parâmetros relativos 
à carga e, por conseguinte, interferindo com a parte de transmissão do sistema. 
Considerando que para haver sobretensão o pólo da função de transferência 

e Q(S> 


-fd(S) 


deve ser negativo, têm-se as seguintes hipóteses. 


Primeira hipótese, com base na equação 3.137: 

1 - <0 e 1 - cü 2 L ’(jC> 0 

p 1 

1-üü^L h C<0 :. - <wL h 

d (oC d 

1-cü 2 L' h C>0 :. —>(jüL' h 
d üjC d 


Ou seja 
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Segunda hipótese, com base na equação 3.137: 


1 - u) 2 L d C > O 

1-u 2 L' d C<0 
Ou seja 


1 - u) 2 L d C > 0 e 1 - u) 2 L ' d C < 0 


—— > 6l»L w 

Co C d 


< 6 oL' h 
ujC d 


LüL d < ~t^C < UjL 'd (140) 

Como já é sabido, em termos de máquina síncrona, a impedância de regime 
permanente ou síncrona Xd será sempre maior que a impedância a essa 
associada de regime transitório X'd, fato que inviabiliza a segunda hipótese. 
Estendendo a análise para o caso em que se foque o eixo d a seguinte 
condição deverá ser verificada: 


1 , 

iüL'< -< (ú L„ 

9 6 oC d 

A equação 3.139 estipula a faixa de valores para a reatância capacitiva 
relativamente aos valores das reatâncias de eixo direto transitória e de regime 
permanente onde a possibilidade de ocorrência do fenômeno de autoexcitação 
no enlace do fluxo de eixo direto é passível de ocorrer. 

De modo semelhante, a equação 3.140 estipula a faixa de valores para a 
reatância capacitiva relativamente aos valores das reatâncias de eixo 
quadratura transitória e de regime permanente para o caso em questão. 

O fenômeno vem então associado ao crescimento exponencial da tensão de 
excitação da máquina síncrona o que pode trazer danos severos ao conjunto 
elétrico independente da atuação rápida do sistema automático de excitação. 



3.3. Considerações Finais 


Pelo que foi exposto, torna-se importante o conhecimento dos parâmetros que 
vão representar a carga que estará conectada a máquina síncrona. Existem 
situações em que determinados arranjos são possíveis em função da abertura 
de determinadas linhas de transmissão ou mesmo situações de baixa carga em 
que a carga capacitiva torna-se relevante. Nesses casos é interessante 
analisar todos os arranjos possíveis determinando o valor da carga capacitiva 
equivalente no sentido de se evitar o fenômeno da autoexcitação. Outro ponto 
que também deverá ser observado são os esquemas em que são utilizados 
compensadores capacitivos no sentido de controlar os níveis de tensão. Tais 
compensadores, dependendo do tamanho, podem concorrer diretamente para 
culminar com o fenômeno da autoexcitação. 
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